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EXECUTIVE  SUMMARY 


Subsurface  contamination  by  organic  compounds  has  become  a  major  environ¬ 
mental  problem  at  various  U.S.Army  installations.  The  technology  exists  to  experi¬ 
ment  with  biological  in-situ  degradation  methods  of  treatment  and  reclamation,  but 
the  factors  controlling  the  chemical  and  physical  environments  of  biodegradation  are 
not  well  known  and  few  site-specific  case  studies  are  documented  in  the  engineering 
literature.  The  purpose  of  this  study  was  to  identify  applicable  in-situ  groundwater 
treatment  methods  and  performance  criteria,  and  to  help  evaluate  these  techniques, 
especially  with  respect  to  limitations  imposed  by  hydrogeologic  conditions. 


Several  in-situ  biorestoration  scenarios  are  described  in  the  report,  mostly  with 
application  to  contamination  in  extensive  sand  aquifers.  Evaluation  of  the  various 
methods,  however,  is  likely  to  be  site  specific,  given  the  complexities  of  most  geologic 
settings  and  groundwater  systems.  The  application  of  in-situ  schemes  for  clayey 
deposits  or  fractured  rock  are  unknown.  A  major  finding  from  this  research  is  the 
potential  effect  of  biofouiing  on  the  site  hydrology.  Permeability,  for  example,  can  be 
reduced  by  several  orders  of  magnitude  due  to  microbial  growth  in  granular  porous 
media.  Laboratory  benchtop  experiments  confirmed  the  permeability  reduction  at 
a  macroscopic  scale.  Mathematical  models  are  presented  also  to  demonstrate  the 
possible  effect  at  the  aquifer  scale,  but  only  site-specific  pilot  studies  can  confirm 
these  results.  .  '  »ir 
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1.  Introduction 


Subsurface  contamaination  of  groundwater  by  organic  compounds  has  become 
a  major  environmental  problem  at  various  U.S.  Army  installations  throughout  the 
country.  It  is  essential  to  find  appropriate  reclamation  methods  to  limit  and  re¬ 
duce  the  extent  of  contamination.  Restoration  of  polluted  aquifers  by  means  of 
biological  in-situ  methods,  which  use  microorganisms  to  degrade  toxic  organic  con¬ 
taminants,  has  recently  been  applied  with  some  success  and  has  been  shown  to  be 
cost-effective  (Wilson  et  al.,  1985,  Wagner  et  al.,1986).  These  in-situ  methods  will 
herein  be  referred  to  in  this  text  as  bioreclamation  or  biorestoration.  The  basic 
method  of  bioreclamation  is  to  enhance  bacterial  growth  in  the  subsurface  by  the 
injection  of  bionutrients.  Heterotrophic  bacteria  can  use  toxic  organic  compounds 
as  a  source  of  carbon  as  well  as  a  source  of  energy,  thereby  degrading  and  detoxi¬ 
fying  the  contaminants.  The  mechanism  of  these  biorestoration  methods,  however, 
is  still  under  research  and  little  is  known  about  the  biochemical  interactions  during 
biodegradation  processes  (Finder  and  Gray,  1984,  Wilson  et  al.,  1985). 

A  further  unknown  factor,  which  will  be  critical  in  actual  field  installations, 
is  the  physical  effect  of  microbial  growth  on  the  hydraulic  properties  of  porous 
media.  The  permeability  of  a  sand  aquifer,  for  example,  can  be  altered  by  orders 
of  magnitude  due  to  microbial  growth.  This  type  of  clogging  effect  would  have  a 
substantial  influence  on  the  transport  of  contaminants  in  groundwater  and  therefore 
on  the  design  of  restoration  schemes. 

2.  Objectives  of  Research 

In  this  research  we  have  tried  to  quantify  the  changes  in  hydraulic  parameters 
due  to  microbial  activity  and  to  outline  the  importance  of  these  changes  for  the 
reclamation  design  and  evaluation  of  performance  criteria.  Our  method  of  analysis 
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included: 


•  A  literature  review  to  identify  possible  in-situ  reclamation  schemes. 

•  Collection  of  data  from  laboratory  permcameter  experiments  to  determine 
changes  in  permeability  under  a  variety  of  controlled  conditions. 

•  Development  of  a  hydrogeologic  simulation  model  capable  of  mathematically 
modeling  two-dimensional  mass  transport  in  groundwater  flow  and  incorporat¬ 
ing  the  effects  of  biological  processes  on  permeability  to  assess  the  importance 
of  these  effects  for  bioreclamation  designs. 

The  goal  of  this  research  will  be  eventually  to  apply  the  experimental  data  and 
the  numerical  model  to  hydrogeologic  settings  typical  of  site  conditions  at  U.S.Army 
installations  for  the  development  of  simulations  of  aquifer  restoration  in  order  to 
evaluate  the  reclamation  procedure  and  its  effects  on  hydraulics  of  the  groundwater 
system. 

3.  Bioreclamation 

IVeatment  Methods 

There  are  a  wide  range  of  possible  treatment  methods  available  to  reclamate 
contaminated  groundwater.  The  conventional  methods  are: 

•  elimination  of  the  contamin2mt  source  by  excavating  it  and  depositing  it  at 
hazardous  waste  disposal  sites, 

•  isolating  the  contaminant  source  by  surrounding  It  with  low  permeability 
slurry  walls,  or  through  the  development  of  groundwater  divides, 
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•  removal  of  contaminated  water  through  pumping  followed  by  above  ground 
treatment, 


Excavating  the  contaminant  source  is  not  always  an  ecomunically  feasible  method. 
The  excavation  costs  are  high  and  workers  aind  the  environment  have  to  be  protected 
against  the  hazardous  chemicals  being  exposed  during  excavation.  The  problem 
sometimes  is  not  resolved  by  depositing  the  contaminant  at  a  hazardous  waste 
disposal  site.  The  construction  of  low  permeability  slurry  walls  can  be  cost-extensive 
(Zipfel  and  Geldner,  1985),  but  isolating  the  contaminant  does  not  eliminate  it. 
There  exists  also  the  possibility  that  certain  chemicals  could  attack  the  slurry  wails 
and  render  them  more  permeable.  Thus,  a  continuous  monitoring  system  becomes 
neccessary.  Removal  of  contaminated  water  through  weU  pumping  is  also  cost- 
effective  for  small  spills,  but  this  becomes  impractical  or  too  expensive  for  large 
contaminant  plumes  in  aquifers. 

In-situ  hydraulic  methods  are  often  used  in  connection  with  chemical  or  biologi¬ 
cal  in-situ  techniques.  If  they  are  applicable  and  prove  to  be  successful,  these  in-situ 
techniques  would  have  the  advantage  to  eliminate  the  contamination  at  place.  The 
site  would  need  to  be  monitored  only  for  the  limited  time  of  treatment.  The  haz¬ 
ardous  chemical  is  left  in  the  ground  and  does  therefore  not  present  a  health  threat 
to  workers.  The  aquifer  can  be  used  again  after  successful  reclamation  (Nagel  et 
al.,  1982).  Among  the  in-situ  methods  bioreclamation  methods  seem  to  show  the 
most  promise  (Wagner  et  al.,  1986). 

Bioreclamation 

The  principle  of  bioreclamation  is  based  on  the  biodegradabiiity  of  the  organic 
contaminant.  Bacterial  growth  is  enhanced  at  the  location  of  the  contaminants  by 
injecting  nutrients  and  oxygen.  Figure  1  shows  a  possible  scenario  of  a  bioreclama¬ 
tion.  The  nutrient  water  is  either  injected  by  a  well  or  infiltrated.  Oxygen  can  be 

3 


! 


Figure  1:  Conceptual  model  for  in  situ  biorratoration.  Nutrients  can  either  be 
injected  by  wells  or  infiltrated  through  flooding  or  spray  irrigation.  Oxygen 
can  either  be  added  in  line  to  the  nutrient  solution  or  directly  injected  into 
the  contaminated  area  through  aeration  wells,  (from  Wagner  et  al.,  1986) 
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added  in-line  or  directly  injected  by  aeration  wells  into  the  contamin2mt  plume. 

Conditions  for  Biodegradation 

Numerous  researchers  have  shown  in  laboratory  experiments  (Alexander,  1981; 
Kobayashi,  1982;  Bouwer,  1984)  that  hazardous  organic  contaminants  are  actually 
transformed  by  bacteria  into  non-toxic  compounds.  Most  compounds  we  degraded 
aerobically  whereas  tower-molecular  weight  halogenated-hydrocarbons  degrade  only 
anaerobically  (Bouwer,  1981;  Wagner  et  al.,  1986).  The  potential  for  utilizing 
biodegradation  in  the  subsurface  also  exists  because  of  recent  work  showing  the 
presence  of  bacteria  (Sykes  et  at.,  1982;  Schwarzenbach  et  al.,  1983;  Wilson  et  al., 
1983).  Bacteria  adapted  to  the  contaminant  in  the  laboratory  can  be  introduced 
into  the  aquifer  but  in  most  bioreclamation  cases  bacteria  indigenous  to  the  aquifer 
are  used. 

Nutrients  must  be  added  to  the  aquifer  in  order  to  enhance  bacterial  activity. 
The  composition  of  these  nutrients  can  be  tailored  for  existing  geochemical  and 
biological  conditions.  In  some  cases,  when  the  contaminant  is  either  toxic  to  bacteria 
or  its  concentration  is  too  low,  bacterial  growth  cannot  be  supported  by  the  organic 
contaminant  alone.  Another  carbon  source  such  as  glucose  or  citrate  has  to  be 
provided  as  a  primary  substrate  to  encourage  bacterial  growth  (Wagner  et  al.,  1986). 

Sufficient  oxygen  supply  is  crucial  for  aerobic  degradation,  where  oxygen  serves 
as  electron  acceptor.  The  injection  of  air  into  the  groundwater  has  limited  applica¬ 
tion  for  contaminations  of  low  concentration,  beca\ise  of  the  relative  low  solubility 
of  oxygen  in  water.  This  difficulty  is  overcome  by  the  use  of  pure  oxygen  and  hydro¬ 
gen  peroxide.  Hydrogen  peroxide  has  the  additional  advantage  of  reducing  clogging 
around  injection  wells  due  to  its  relative  toxicity  to  bacteria,  and  it  decomposes 
quickly  into  water  and  oxygen  thus  reducing  its  toxicity  and  making  oxygen  avail¬ 
able  to  the  bacteria. 
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Design  of  Bioreclamation  Systems 

Nutrients  and  oxygen  have  to  be  injected  or  transported  into  the  contamination 
site.  Injection  and  extraction  systems  are  the  two  basic  units  of  a  biorestoration 
scheme.  The  transport  of  nutrients  and  oxygen  can  be  gravity  driven  or  forced  by 
pumping  wells.  If  the  soil  is  permeable  enough  to  allow  percolation  the  nutrient 
water  can  be  applied  to  the  contamination  site  by  spray  irrigation  or  flooding.  In 
less  permeable  soils  subsurface  drains  might  be  appropriate  (Figure  2).  The  ap¬ 
plication  of  subsurface  drains  is,  however,  limited  to  shallow  contaminant  plumes 
in  unconsolidated  media  because  the  installation  of  the  drainage  system  to  greater 
depths  could  be  too  cost-intensive  (Wagner,  1986).  A  wide  variety  of  arrangements 
of  injection  and  extraction  wells  are  possible.  In  general,  rows  of  wells  injecting  nu¬ 
trients  and  oxygen  are  arranged  upstream  of  the  contamination  plume.  Extraction 
wells  are  located  downstream  of  the  contamination  site  to  collect  the  water.  The 
water  can  be  further  treated  aboveground  and  recharged  into  the  aquifer  through 
the  injection  wells  after  addition  of  nutrients  and/or  oxygen.  Another  group  of 
wells,  which  inject  fresh  groundwater  and  are  located  around  the  wells  injecting 
the  recharged  water  with  nutrients,  create  a  hydraulic  mound  that  prevents  further 
contamination  of  the  aquifer  by  recharge  water  (Figure  3).  Additional  recharge 
wells  arranged  along  the  contamination  plume  could  optimize  nutrient  and  oxygen 
supply  to  the  contaminated  area  (Figure  4).  If  the  groundwater  flow  direction  and 
the  extent  of  the  plume  is  well  known,  a  row  of  injection  wells  perpendicular  to  the 
groundwater  flow  direction  infiltrating  nutrients  and  oxygen,  can  eventually  inter¬ 
cept  the  spread  of  the  contamination  (Figure  5).  In  all  cases  monitoring  wells  have 
to  control  the  spatial  extent  of  the  contamination  to  make  sure  that  the  operational 
system  is  working  efficiently  and  that  the  contamination  plume  is  not  spreading 
further. 
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Figure  2;  Bioreclaunation  with  a  subsurface  drain.  In  less  permeable  soils  a  sub¬ 
surface  drain  may  be  more  appropriate  than  extraction  wells,  (from  Wagner 
ei  al.,  1986) 
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Equipoleatikl  Lines 
Limit  of  Recharge  Area 
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Recharge  Water 
Clean  Infiltration  Water 


®  Clean  Water  Infiltration 

•  Recharge  Water  Infiltration 

®  Extraction  Well 

•  Obeerration  Well 


Figure  3:  Bioreclamation  with  fresh  water  infiltration.  During  in-situ  bioreclama¬ 
tion  the  polluted  extracted  groundwater  is  often  treated,  and  after  addition 
of  nutrients  and  oxygen,  recharged  into  the  subsurface.  Wells  around  the 
recharge  wells,  inject  clean  groundwater  to  prevent  a  contamination  of  the 
aquifer  with  the  polluted  recharge  water,  (after  Battermann  and  Werner, 

1984) 
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Flow  Direction 


o  Extraction  Well 

*  Recharge  Water  Infiltration 


Figure  4:  Bioreclamation  with  additional  injection  wells.  This  figu’e  shows  a 
possible  configuration  of  injection  and  extraction  wells  suggested  by  Valocchi 
and  Rittmann  (1986).  Injection  wells  along  the  contaminated  plume  would 
facilitate  transport  of  nutrients  and  enhance  biodegradation. 
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Figure  5;  Bioreciamation  with  injection  wells  downstream  from  the  contaminant 
source.  A  possible  configuration  scheme  for  a  bioreciamation  is  to  Intercept 
the  spread  of  a  contaminant  plume  through  injection  wells  downstream  from 
the  contaminant  source,  (after  Bouwer,  1984) 


10 


4.  Performance  Criteria 


The  following  performejice  criteria  for  the  successful  implementation  of  a 
bioreclamation  were  identified: 

•  previous  experiences 

•  operation  and  monitoring  requirements 

•  costs 

•  biodegradability 

•  bacterial  contamination 

•  hydrogeological  implications 

Experiences  from  previous  bioreclamation  applications  show  that  many  organic 
contaminants  can  be  biodegraded  to  concentrations  below  detection  level  (Nagel, 
1982;  Jhaveri,  1983).  Bioreclamation  has  been  applied  only  to  spill  sites  so  far, 
but  since  there  is  evidence  of  adapted  bacterial  communities  below  waste  disposal 
sites,  Wagner  et  al.(l986)  suggest  that  bioreclamation  methods  are  applicable  to 
such  cases  as  well.  Depending  on  the  duration  of  the  bioreclamation  system,  opera¬ 
tion  and  maintenance  requirements  can  be  quite  expensive.  A  series  of  biochemical 
parameters  such  as  pH,  dissolved  oxygen,  concentration  of  contaminant,  bacterial 
numbers,  etc.,  have  to  be  monitored  constantly  in  order  to  maintain  the  efficiency 
of  the  reclamation  process.  These  puameters  can  be  corrected  by  addition  of  ap¬ 
propriate  nutrients,  oxygen,  or  by  incubation  of  bacteria.  In  many  bioreclamation 
cases,  clogged  wells  had  to  be  treated.  The  use  of  hydrogen  peroxide  has  shown  to 
be  an  effective  agent  for  treating  clogging  effects  (Yaniga,  1984).  Estimations  for  op¬ 
erational  and  permanent  costs  are  difficult  to  generalize  because  of  the  site-specific 
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characteristics  of  a  contamination  site.  It  is  also  difficult  to  make  any  generalizations 
about  costs  since  the  success  of  a  bioreclamation  method  is  dependant  on  so  many 
different  paraimeters.  Wagner  et  a!.(1986)  gives  one  example  of  a  contamination 
site  where  the  in-situ  treatment  of  the  contaminated  groundwater  through  biorecla¬ 
mation  was  an  order  of  magnitude  less  costly  than  the  above  ground  treatment  of 
the  polluted  groundwater.  Laboratory  experiments  have  to  prove  that  the  organic 
contaminant  is  biodegradable  by  the  bacterial  population  present.  In  some  cases, 
a  high  heavy  metal  concentration  could  be  toxic  to  bacteria  aoid  prevent  successful 
biodegradation.  The  possibility  also  exist  for  contaminants  to  be  transformed  by 
bacteria  into  still  more  toxic  compounds  (Alexander,  1981).  The  rate  of  degrada¬ 
tion  depends  on  the  characteristics  of  the  organic  contaminant  and  the  bacterial 
population.  If  the  degradation  requires  a  long  time,  say  5-6  years,  to  occur,  then 
due  to  high  operational  costs  the  expenses  for  bioreclamation  might  exceed  those  of 
a  more  conventional  water  treatment  method.  In  other  cases  there  exists  concern 
about  bacterial  contamination.  Little  is  known  about  the  transport  of  bacteria  in 
aquifers.  Some  tracer  studies  with  viruses  indicate  that  microorganisms  can  travel 
hundreds  of  meters  (McDowell-Boyer  et  al.,  1986).  The  proximity  of  water  supply 
wells  could  raise  concern  about  bioreclamation  methods  where  bacterial  growth  is 
enhanced.  In  addition  to  these  restrictions  the  success  of  a  bioreclamation  scheme 
will  depend  very  strongly  on  the  hydrogeological  situation  of  an  aquifer,  as  described 
in  the  next  section. 

S.  Hydrogeological  Implications 

Most  of  the  known  in-situ  bioreclamation  cases  have  b'ien  applied  to  highly  per¬ 
meable  aquifers.  Contamination,  however,  can  occur  in  all  types  of  media,  ranging 
from  soils  and  clays  to  gravel  and  fractured  rock.  The  tranport  of  nutrients  to  the 
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contaminated  groundwater  zones  is  dependent  on  the  fluid  transport  a.id  the  reac¬ 
tion  rate  of  the  nutrient.  In  geologic  media  of  low  permeability  the  fluid  velocities 
are  very  low.  If  the  reaction  rate  is  fast  relative  to  the  fluid  velocity,  the  nutrient 
is  likely  to  be  used  up  by  bacteria  closer  to  the  nutrient  source.  Under  these  con¬ 
ditions  the  nutrient  would  not  be  transported  eflectively  to  a  large  contaminated 
area  and  little  or  no  biodegradation  may  take  place.  Biological  clogging  around  the 
nutrient  source  is  very  likely  to  occur.  If  the  reaction  rate  was  slow  compared  to  the 
advection  of  solute  then  it  would  be  possible  that  the  nutrient  would  be  transported 
to  the  contaminated  area,  even  in  low  permeability  media.  The  ratio  of  advective 
transport  and  reactive  transport  b  reflected  in  the  dimensionless  Damkoehler  num¬ 
ber  which  could  serve  as  an  indicator  of  the  feasibility  of  a  bioreclamation  method 
in  an  media  of  certain  permeability.  The  distance  from  the  nutrient  source  to  the 
contamination  site  in  lower  permeability  media  has  to  be  reduced  as  much  as  pos¬ 
sible.  Injection  well  spacings  should  therefore  be  reduced  and  the  use  of  subsurface 
drains  incorporated  if  appropriate.  Another  concern  in  low-permeability  aquifers  is 
the  reduction  of  permeability  due  to  clay  dispersion.  Wagner  et  al.  (1986)  reports 
of  soil  conditioners  which  supposedly  prevent  clay  dispersion  and  therefore  enhance 
therefore  transport  of  nutrients  and  bacteria. 

The  groundwater  velocity  and  flow  direction  can  be  imposed  in  general  by 
the  arrangement  of  the  recharge/discharge  welb  and  their  pumping  rates.  Dur¬ 
ing  bioreclamation  cases,  however,  this  imposed  flow  pattern  b  subject  to  changes 
due  to  permeability  reductions.  The  addition  of  nutrients  and  oxygen  source  in  the 
recharge  water  b  responsible  for  permeability  reductions  for  numerous  reasons; 

•  If  the  iron  or  manganese  content  in  soils  or  sediments  b  high,  precipitation  of 
iron  and  manganese  results  from  oxidation  in  the  aquifer  (Smith  and  Tuovinen, 
1985;  Hackett,  1987). 
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•  A  sediment  with  a  high  clay  content  cam  bs  subject  to  clay  dispersion,  amd 
therefore  to  permeability  reductions  if  leaurhed  with  a  sodium-rich  nutrient 
solution 

•  The  various  injection  methods  for  supplying  oxygen  to  the  subsurface  involve 
the  possibility  that  oxygen  in  the  form  of  bubbles  is  created  inducing  a  multi¬ 
phase  flow  regime.  If  hydrogen  peroxide,  for  instamce,  transforms  very  rapidly, 
oxygen  may  bubble  out  as  a  separate  fluid  phase. 

•  Bau:terial  growth  is  enhanced  by  the  addition  of  nutrients,  and  bacterial  biomass 
might  increase  extensively  thus  reducing  permeability. 

Depending  on  the  significance  of  the  above  mentioned  effects,  permeability  re¬ 
ductions  in  biochemically  active  zones  can  alter  flow  velocity  and  Sow  direction 
significantly.  Through  the  induced  heterogeneities  due  to  the  created  low  permeabil¬ 
ity  zones  the  dispersion  of  the  contaminant  might  also  increase.  As  a  consequence, 
the  prediction  of  groundwater  flow  and  contaminant  transport  mav  be  miscalcu¬ 
lated  if  the  above  mentioned  effects  are  not  considered.  Despite  the  fact  that  these 
phenomena  are  difficult  to  quantify,  their  consideration  for  successful  planning  of 
bioreclamation  methods  seems  mandatory. 

6.  Literature  Review 

Many  laboratory  experiments  have  employed  column  permeameters  to  observe 
changes  in  hydraulic  conductivity  or  permeability  reaulting  from  microbial  activ¬ 
ity  (Allison,  1947,  McCalla,  1950,  Gupta  and  Swartiendruber,  1962,  Poulovassilis, 
1972,  Frankenberger  et  al.  1979,  Shaw  et  al.,  1985).  Hydraulic  conductivity  is  a 
parameter  used  to  characterize  the  ability  of  a  porous  medium  to  transmit  fluids 
for  a  given  hydraulic  gradient.  It  is  a  volume-sveraged  parameter  defined  as  the 
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proportionality  constant  in  Darcy’s  law; 


q^-K 


ds 


Where:  h  is  the  hydraulic  head,  Im] 

s  is  the  direction  coordinate  of  Sow,  [m] 
q  is  the  specific  discharge,  [m/s] 

K  is  the  hydraulic  conductivity,  [m/s] 

The  hydrrulic  conductivity  is  related  to  the  permeability  k  through: 

kjg 


K  = 


(1) 


(2) 


Where:  p  is  the  fluid  density,  [A:j/m*] 

fi  is  the  fluid  viscosity,  [^s/m*] 
g  is  the  acceleration  of  gravity,  [m/s*] 
k  is  the  intrinsic  permeability,  [m*] 

The  permeability  of  a  porous  medium  is  controlled  by  lithologic  factors  such  as 
grain  packing,  pore  shape,  and  specific  surface,  tortuosity,  and  connectivity  (Bear, 
1972).  For  many  geologic  media  the  permeability  is  also  greatly  affected  by  the 
degree  of  fracturing  and  cementation. 


Figure  6  gives  a  summary  of  the  references  cited  here  which  examine  the  influence 
of  microbial  growth  on  permeability  under  biotic  and  abiotic  conditions.  A  perusal 
of  the  literature  has  shown  that  three  types  of  experimental  approach  have  been 
used  to  inhibit  biological  activity  within  the  columns: 


•  leaching  under  sterile  conditions  (Allison,  1947;  Poulovassilis,  1972;  Franken- 
berger  et  al.,  1979;  Shaw  et  al.,  1985). 

•  leaching  with  toxic  solutions  (Allison,  1947;  Gupta  and  Swartzendruber,  1962; 
Shaw  et  al.,  1985). 
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Figure  6;  Literature  review  of  laboratory  experiments  evaluating  the  influence  of 
bacteria  on  the  hydraulic  conductivity  of  porous  media. 
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•  leaching  under  conditions  with  temperature  slightly  above  freezing  (McCalla, 
1950;  Gupta  and  Swartzendruber,  1962). 

No  changes  in  permeability  could  be  observed  in  experiments  where  biological  ac¬ 
tivity  in  the  porous  media  was  inhibited.  Non-sterilized  samples  ,  however,  showed 
a  substantial  decrease  in  hydraulic  conductivity.  Frankenberger  et  al.  (1979) 
compared  the  effect  of  different  nutrients  on  the  hydraulic  conductivity.  His  re¬ 
sults  show  that  a  solution  of  glucose  and  KNO3  was  the  most  effective  in  terms 
of  hydraulic  conductivity  reduction.  Okubo  and  Matsumoto(1983),  and  Shaw  et 
al.(l985),  showed  that  the  permeabilty  was  still  more  reduced  when  their  leaching 
solution  contained  suspended  solid  particles. 

The  results  of  these  experiments  document  the  ability  of  microbial  processes  to 
reduce  the  permeability  of  a  granular  porous  medium  by  as  much  as  two  or  three 
orders  of  magnitude  ,  at  least  for  periods  of  reactive  fluid  flow  on  the  order  of 
hundreds  of  days.  This  is  clear  evidence  that  microbial  degradation  effects  on  the 
hydraulic  properties  ought  to  be  taken  into  account  in  any  quantitative  assessment 
of  aquifer  restoration  by  in  situ  methods.  One  problem,  however,  is  deciding  on 
the  means  of  quantifying  the  permeability  reductions.  Gupta  and  Swartzendruber 
(1962)  sought  a  relationship  between  hydraulic  conductivity  and  the  total  mass  of 
microorganisms.  From  their  experiments,  changes  in  hydraulic  conductivity  were 
observed  only  above  a  threshold  level  of  4x10*  bacteria  per  gram  of  soil.  Franken¬ 
berger  et  al.(l979)  observed  a  more  uniform  decline  in  permeability  with  increase  in 
bacterial  mass,  which  they  gaged  with  phosphatase  production.  Okubo  and  Matzu- 
moto  (1983)  observed  in  their  experiments  an  increase  of  clogging  when  either  the 
substrate  input  concentration  or  the  suspended  particle  concentration  or  both  of 
them  were  increased.  Furthermore  they  could  show  that  there  exist  three  distinct 
zones  exhibiting  different  clogging  rates.  The  first  rone  was  identified  as  the  aerobic 
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phase  and  showed  the  fastest  growth  rate.  The  second  zone  was  a  transition  zone 
between  the  aerobic  and  anaerobic  phase,  showing  a  slight  increase  in  permeabil¬ 
ity.  Finally  a  third  zone,  the  anaerobic  phase,  was  characterized  by  a  continuous 
decrease  of  permeability,  but  at  a  lower  rate  than  during  aerobic  growth. 

It  is  not  yet  clear  why  or  how  the  bacterial  activity  causeti  the  permeability 
reduction.  The  following  factors  may  all  be  of  importance  : 

•  mass  of  bacteria 

•  solid  products  produced  by  bacteria  like  exocellular  polymers 

•  gases  produced  by  bacteria 

•  disaggregation  of  soil  grains  in  the  porous  media  through  bacterial  consump¬ 
tion  of  aggregating  organic  compounds 

•  inorganic  dissolution/precipitation  reactions  induced  by  bacteria 

Considering  that  the  threshold  level  of  4x10*  bacteria  per  gram  sand  (Gupta  and 
Swartzendruber,  1962)  represents  only  0.0003  %  of  the  total  volume  of  one  gram 
volume  of  sand,  bacteria  mass  may  not  be  the  reason  for  permeability  reduction 
alone.  Micrographs  made  by  Shaw  et  al.(l985)  show  in  their  experiments  extra¬ 
cellular  polymers  as  a  main  cause  for  bacterial  clogging  processes.  McCalla  (1950) 
and  Poulovassilis  (1972)  identified  gas  production  as  a  major  reason  for  permeabil¬ 
ity  decreases.  Other  factors  such  as  inorganic  reactions  (Quirk  and  Schofield,  1955; 
Suarez  et  al.,  1984)  and  swelling  and  shrinking  of  clays  (Goldenberg  et  al.  ,  1983; 
Fernandez  and  Quigley,  1985)  may  also  have  played  a  role  in  affecting  permeability 
in  the  porous  media  used  in  the  experiments. 

Observations  in  the  field  suggest  that  bacteria  are  responsible  for  permeabilty 
reductions  in  the  subsurface.  Sealing  effects  at  infiltration  ponds  during  waste  wa- 
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ter  recharge,  as  well  as  during  river  infiltration  were  reported  by  Geldner  (1980) 
and  Schwarzenbach  et  al.(1983).  Numerous  examples  of  well  clogging  are  also  re¬ 
ported  (Ehrlich,  et  al.,  1973;  Hijnen  and  van  der  Kooij,  1984;  Smith  and  Tuovinen, 
1985).  In  many  cases  these  phenomena  could  be  related  to  bacterial  activity.  Field 
groundwater  samples  taken  by  Harvey  et  al.(l984)  in  a  aquifer  contaminated  by  a 
sewage  treatment  effluent,  indicated  orders  of  magnitude  higher  number  of  bacteria 
within  the  contamination  plume. 

7.  Laboratory  Experiment 

With  the  exception  of  Shaw  et  ai.(l985)  who  performed  their  experiments 
with  sintered  glass  beads,  all  of  the  above  cited  studies  have  used  a  porous  medium 
composed  of  geologic  material  such  as  sand  or  loam.  Because  of  the  inherent  min- 
eralogical  complexity  of  a  geologic  material  it  is  difficult  to  isolate  the  relative 
importance  of  the  numerous  biological,  chemical,  and  physical  processes  that  could 
affect  permeability  changes.  All  that  is  known  so  f2U'  is  that  permeability  reductions 
in  shallow  aquifers  mainly  occur  as  a  result  of  bacterial  activity. 

Attempts  to  quantify  these  effects  with  respect  to  bacterial  mass  were  done  by 
Frankenberger  et  al.(1979).  Frankenberger  used  a  clayey  loam  in  his  experiment 
as  the  porous  medium.  It  is  doubtful,  however,  that  Frankenberger’s  results  can 
be  generalized  for  other  geological  materials.  The  biological  system  of  most  experi¬ 
ments  consisted  of  a  mixed  system  of  aerobic  and  anaerobic  growth.  Bacterial  gases 
formed  under  anaerobic  conditions  during  the  denitrifying  and  methanogenic  phase 
(Figure  7)  might  have  a  major  impact  on  permeability  changfa. 

In  order  to  prove  that  permeability  reductions  are  due  mainly  to  bacterial  gas 
production,  Poulovassilis  degassed  his  permesimeter  column  under  vacuum  and  ob¬ 
served  a  significant  permeability  increase.  Applying  a  vMuum  onto  the  perme- 
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Figure  7:  Possible  gas  production  by  bacteria.  A  microbial  population  is  charac¬ 
terized  by  its  electron  acceptors.  Depending  on  the  electron  acceptor  present 
different  gases  can  be  produced  by  bacteria.  Nitrogen  and  methauie  gM  are 
likely  to  be  formed  under  denitrifying  and  methanogenic  conditions. 
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ameter,  however,  could  have  destroyed  a  biological  matrix  formed  during  bacteria! 
growth  which  could  also  have  exerted  an  influence  on  the  permeability.  It  is  ques¬ 
tionable  therefore  to  what  extent  the  permeability  increase  can  be  ascribed  to  gas 
release,  since  degassing  could  have  disturbed  the  porous  matrix.  The  above  workers 
neither  offered  an  explanation  for  the  hydraulic  conductivity  reduction,  nor  did  they 
give  a  quantification  of  hydraulic  conductivity  reduction  due  to  biological  processes. 

Objectives  of  Laboratory  Column  Experiments 

A  basic  goal  of  this  study  was  to  overcome  deficiencies  of  the  earlier  work  cited 
above,  so  as  to  help  evaluate  the  role  of  permeability  changes  in  in-situ  reclamation. 
We  want  to  establish  with  our  laboratory  experiments  a  quantifiable  relationship 
between  permeability  and  substrate  concentration.  Glass  beads  were  used  as  the 
porous  medium  in  the  experiments,  to  exclude  all  non-biclogical  factors  on  perme¬ 
ability  change,  The  substrate  concentration  and  defined  electron  acceptor  conditions 
are  a  means  of  controlling  biomass.  By  maintaining  a  heterotrophic,  aerobic  biolog¬ 
ical  environment,  the  only  gas  which  could  form  due  to  bacterial  activity  is  likely  to 
be  COj.  The  solubility  of  CO2  under  atmospheric  pressure,  however,  is  sufficiently 
high  that  gas  formation  would  not  be  expected  to  be  significant  with  the  lov  sub¬ 
strate  concentrations  used  in  the  experiments.  Hence,  gas  evolution,  was  excluded 
as  a  permeability  reducer  for  the  experimental  conditions  chosen.  The  permeability 
changes  were  monitored  along  the  length  of  the  column  and  over  time  for  a  constant 
fluid  flow  velocity  and  a  given  substrate  input  concentration.  Measurements  of  the 
substrate  concentration  along  the  column  as  well  as  theoretical  calculations  with  a 
biofilm  model  were  used  to  help  evaluate  the  effects  of  transport  on  substrate  con¬ 
centration  in  a  flow  field.  The  substrate  distribution  is  related  to  bacterial  mass  by 
means  of  a  biofilm  model.  Microscopic  analysis  after  each  ixperiment  confirmed  the 
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theoretical  predictions  by  the  biofilm  model  concerning  biomass  distribution  made. 
Hence  permeability  changes  can  be  directly  related  to  the  distribution  of  substrate 
and  bacterial  mass  within  the  column.  Permeability  distributions  are  compared  for 
different  substrate  input  concentrations. 

Microscopic  examination  of  the  biomass  formed  can  give  some  evidence  for  the 
cause  of  permeability  chemges.  It  provides  some  insight  into  the  shape  and  size 
of  the  biomass.  This  may  establish  if  the  decrease  in  hydraulic  conductivity  is 
primarily  due  to  volume  reduction  of  the  void  space  or  due  to  the  alteration  of 
surface  characteristics  of  the  solid  matrix. 

Another  hydraulic  parameter  which  is  addressed  in  our  research  is  the  hydrody¬ 
namic  dispersion  coefficient.  Values  of  the  dispersion  coefficient  measured  in  lab¬ 
oratories  generally  differ  greatly  form  those  evaluated  in  field  studies  (Freeze  and 
Cherry,  1979).  Dispersivity  values,  measured  in  the  laboratory  columns,  however, 
could  provide  at  least  qualitative  information  on  the  influence  of  bacterial  growth 
on  the  dispersion  coefficient. 

Description  of  Laboratory  Apparatus 

A  fixed  wall  permeameter  was  used  to  measure  permeability  (Figure  8).  A 
peristaltic  pump  provides  a  constant  flow  while  manometers  ports  along  the  per- 
mcameter  measure  the  changes  in  hydraulic  head. 

Glass  beads  were  used  as  the  porous  material  because  the  possibility  of  grain 
rearrangement  during  flow  and  adsorption/desorption  effects  are  minimized.  Fur¬ 
thermore,  no  shrinking  or  swelling  as  with  clayey  material  amd  no  mineral  dissolu¬ 
tion  needs  to  be  considered  with  glass  beads.  The  organic  carbon  content  can  be 
controlled  and  the  shape  and  surface  area  of  the  grains  are  well  defined.  A  bead 
diameter  which  is  close  to  a  medium  sand  was  chosen  for  the  first  set  of  experi- 
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Figure  8:  Laboratory  apparatus  for  permeability  experiment.  Nutrients  are  sup¬ 
plied  from  the  supply  bottle  through  the  pump  to  the  permeameter.  The 
permeameter  was  inoculated  initially  with  a  sewage-derived  bacteria  solution. 
Manometers  along  the  permeameter  track  pressure  differences  due  to  perme¬ 
ability  changes. 
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ments.  Glass  beads  of  the  FERRO  Corporation  are  within  a  range  of  0.177  and 
0.210  mm.  No  more  than.  2  %  are  irregularly  shaped  and  90  %  of  them  are  true 
spheres.  The  bulk  specific  gravity  averages  close  to  1.5  and  the  density  of  the  solid 
glass  is  approximately  2500  kg/m^. 

A  photograph  of  the  laboratory  set-up  is  shown  in  Figure  9.  Plexiglass  was  used 
to  fabricate  the  permeameters,  because  the  experimental  solution  does  not  contain 
highly  adsorbable  organic  compounds.  Nylon  fittings  came  from  the  Swagelok  Cor¬ 
poration  .  The  column  length  was  15  cm  and  its  cross-sectional  area  was  15.55 
cm*.  Manometer  or  ;lets  were  placed  along  the  column.  Close  to  the  inflow  they 
were  arranged  in  shorter  intervals,  since  the  highest  permeability  reductions  were 
expected  there.  The  manometer  fittings  intrude  slightly  into  the  porous  medium 
to  avoid  wall  effects.  The  glass  beads  were  placed  on  a  nylon  mesh  with  a  mesh 
opening  of  0.125  mm  (Small  Parts  Inc.).  The  transparent  plexiglass  column  was 
wrapped  with  aluminum  foil  to  prevent  growth  of  photosynthetic  microorganisms. 
It  was  easily  removed  to  observe  possible  changes  of  color  in  the  porous  sample  due 
to  biofilm  growth.  A  multichannel,  mini-c2U'tridge  pump  (ISMATEC,  Model-No. 
NK7624-02)  maintains  constant  column  flow.  This  pump  can  sustain  a  flowrate 
range  of  8.3x10"’°  to  4.4xlO"*m*/s.  Silicon  is  an  inert  material  but  permeable  to 
gases,  and  as  such  it  was  used  for  the  tubing  from  the  supply  bottle  to  the  inlet  of 
the  permeameter  where  oxygen  saturation  is  desired.  Teflon  tubing  which  is  very 
inert  and  less  permeable  to  oxygen  was  used  between  the  permeameter  outflow  and 
sample  bottle. 

Experimental  Procedure 

The  experimental  procedure  is  outlined  in  Figure  10.  The  whole  set-up  was 
sterilized  before  the  experiments  started.  The  glassware  was  sterilized  by  heating 
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Figure  10:  Flow  chart  of  experimental  procedure 


in  an  autoclave.  The  plexiglass  permeameter  and  the  plastic  tubing  were  disinfected 
with  a  10  %  chlorox  solution.  The  column  influent  was  prepared  with  distilled  and 
sterilized  water.  The  gloss  beads  were  poured  into  the  column  filled  with  water 
in  order  to  avoid  air  ii^clusions.  The  column  was  closed  and  the  virgin  hydraulic 
conductivity  measurement  was  madewhich  served  as  reference  value  for  the  initial 
state. 

A  bromide  ion  tracer  was  used  to  provide  an  estimate  of  the  hydrodynamic 
dispersion  coefficient.  The  bromide  breakthrough  response  was  analyzed  with  an 
advection/dispersion  transport  model  (Javandel  et  a].,  1984)  to  evaluate  the  disper¬ 
sion  coefficient.  The  bromide  ion  was  analyzed  by  measuring  the  ion  concentration 
with  an  ion-chromatograph.  The  experimental  break-through  curve  was  fitted  with 
an  analytical  solution  from  which  the  dispersion  coefficient  was  evaluated. 

A  bacteria  solution  was  prepared  by  diluting  primary  settled  sewage.  This  so¬ 
lution  was  pumped  into  the  column.  The  pumping  was  interrupted  for  six  hours  so 
that  bacteria  could  settle  in  the  porous  medium.  A  phosphate  buffer  solution  was 
added  to  the  bacterial  solution  to  prevent  any  cell-lysis. 

After  inoculation,  other  neccessary  nutrients,  such  as  a  phosphate  buffer,  mag¬ 
nesium  sulfate,  calcium  chloride,  ferric  cloride  {Standard  Methods,  1976)  were  added 
and  acetate  was  used  as  the  carbon  source.  Hydraulic  conductivity  measurements 
were  taken  once  or  twice  a  day  in  the  beginning,  and  every  second  day  at  a  later 
stage  of  the  experiment.  The  oxygen  content  of  the  outflow  and  the  pH  was  mea¬ 
sured  once  a  week,  to  ensure  that  the  system  is  aerobic  and  that  no  major  biological 
changes  occurred. 

One  tracer  test  with  radioactive  labeled  acetate  as  the  carbon  source  was  per¬ 
formed  to  evaluate  the  uptake  of  acetate  by  bacteria.  One  liter  of  nutrient  solution 
was  autoclaved  and  20^1  of  radioactive  carbon-14  [acetic  acid-2-^*C}  was  added. 
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so  that  the  reactivity  was  at  about  5000  DPM/ml  (DPM  =  dbintegrations  per 
minute).  The  permeameter  was  leached  with  this  solution.  Samples  were  taken 
from  the  inflow,  outflow  and  along  the  colunm.  The  reactivity  of  the  samples  in¬ 
dicating  the  percentage  of  uptake  of  the  carbon  source  by  bacteria  was  measured 
with  a  liquid  scintillation  analyzer. 

At  the  end  of  the  experiment  a  final  bromide  tracer  test  detected  changes  of 
dispersivity  after  bacterial  growth.  Samples  of  glass  beads  were  taken  along  the 
column  after  the  experiment  and  analyzed  under  a  light  microscope  to  get  visual 
insight  into  volume  and  shape  of  the  bacterial  mass. 

Experimental  Results 

The  hydraulic  conductivity  was  observed  over  time  for  two  input  concentrations 
of  acetate:  1  mgjl  and  '20  mg  (I  (Figure  11).  The  hydraulic  conductivity  was  eval¬ 
uated  after  Darcy’s  law  for  a  column  length  of  10  cm  by  measuring  the  hydraulic 
head  differences  at  the  in  and  outflow  of  the  column  maintaining  a  constant  flux. 
There  exists  for  both  cases  a  lag-phase  during  which  the  hydaulic  conductivity  de¬ 
creased  very  slowly  or  not  at  all,  because  bacterial  growth  had  to  develop  before  it 
showed  a  pronounced  effect  on  hydraulic  conductivity.  A  phase  of  exponential  de¬ 
crease  of  hydraulic  conductivity  followed  the  lag-phase.  The  hydraulic  conductivity 
decreased  about  twice  as  fast  for  the  higher  concentration  case  as  for  the  lower  con¬ 
centration  case.  After  800  hours  the  hydraulic  conductivity  seemed  to  change  little 
for  the  lower  concentration  case.  The  observed  slight  fluctuations  resulted  probably 
from  sloughing  and  replugging  of  biofilm  mass.  For  the  higher  concentration  case 
the  hydraulic  conductivity  kept  decreasing  after  400  hours.  In  both  cases  the  total 
hydraulic  conductivity  decreased  significantly:  in  the  lower  concentration  case  by 
about  two  orders  of  magnitude,  and  in  the  higher  concentration  case  by  about  five 
orders  of  magnitude.  Almost  all  changes  in  hydraulic  head  occurred  within  the 
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Figure  11:  Hydraulic  conductivities  measured  over  time.  They  we  average  val¬ 
ues  for  a  column  of  10  cm  length,  represented  as  the  ratio  of  the  measured 
hydraulic  conductivity  at  time  t  to  its  initial  value  at  time  zero. 
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first  1  cm  behind  the  column  inflow  (Figure  12).  In  the  lower  concentration  case, 
no  change  of  hydraulic  head  above  1  cm  was  observed  at  all.  The  fact  that  most 
of  the  bacterial  activity  occurs  only  in  the  first  1  cm  was  confirmed  by  analyzing 
the  results  of  the  labeled  acetate  measurements  which  were  taken  along  the  column 
during  the  lower  concentration  experiment.  It  shows  that  80  %  of  the  acetate  was 
consumed  within  the  first  1  cm  by  bacteria  (Figure  13). 

Over  the  total  column  length,  99  %  of  the  acetate  was  taken  up  by  bacteria  and 
transformed  into  biomass  and  carbon  dioxide.  The  specific  yield  coeificient  or  cell 
yield  Y  which  represents  the  ratio  of  mass  of  acetate  transformed  into  cell  mass  to 
the  mass  of  acetate  consumed, 

y  _  grams  of  acetate  transformed  into  cell  mass 

grams  of  acetate  consumed  '  ' 

was  evaluated.  The  difference  between  the  mass  of  £w:etate  consumed  and  the  mass 

of  acetate  transformed  into  CO]  corresponds  to  the  mass  of  acetate  transformed 

into  biomass. 

An  average  dissolved  oxygen  concentration  of  about  8  mg /I  at  the  effluent  and 
a  constant  pH  of  7.1  measured  during  the  experiment  of  lower  input  concentration 
showed  that  the  biological  system  is  governed  by  aerobic  heterotrophic  metabolism. 

The  dispersion  coefficient  seems  to  have  increased  slightly  (two-fold)  after  bacte¬ 
rial  growth,  as  shown  in  the  results  of  the  bromide  tracer  test  performed  before  and 
after  the  higher  concentration  experiment  (Figure  14).  This  increase  is  probably 
due  to  the  heterogeneity  and  therefore  increased  tortuosity  of  the  porous  medium 
induced  by  the  bacteria. 

Based  on  these  preliminary  experimental  results,  it  has  become  clear  that  future 
experiments  should  confirm  the  observed  trends: 

•  a  steady  state  hydraulic  conductivity  after  an  exponential  decrease  due  to 
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Figure  12:  Hydraulic  heads  represented  along  the  length  of  the  column  for  different 
times. 
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Figure  14:  Breakthrough  curve  of  a  bromide  tracer  test,  a)  Breakthrough  curve  of  a 
bromide  tracer,  measured  with  no  bacterial  growth  present,  b)  Breakthrough 
curve  of  a  bromide  tracer  with  bacterial  growth  after  400  hours  of  leaching 
with  a  substrate  concentration  of  20  mg/1 
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bacterial  growth 


•  an  empirical  relationship  between  biomass  expressed  in  terms  of  substrate 
concentration  and  hydraulic  conductivity 

•  the  hydraulic  conductivity  reduction  occurs  very  close  to  the  nutrient  source 

•  an  increase  of  dispersivity  due  to  bacterial  growth 

Other  factors  such  as  varying  flow  rates  and  grain  size  of  the  porous  medium  may 
also  affect  the  nutrient  distribution  and  hence  the  biomass  and  hydraulic  conduc¬ 
tivity.  Time  did  not  permit  the  experimental  evaluation  of  these  parameters,  but 
a  theoretical  analysis  of  the  infiuence  of  these  factors  was  performed  with  a  biofilm 
model  . 

8.  Biofilm  Modeling 

The  mathematical  biofilm  model  used  here  to  simulate  the  experimental  set-up 
was  deve!oped  by  Rittmann  (1979).  This  model  calculates  substrate  concentrations 
present  at  any  distance  2und  time  within  a  column  for  a  given  input  concentration  . 
The  biofilm  model  couples  the  one  dimensional  advection-dispersion  equation  with 
a  biological  reaction  term.  The  reaction  term  includes  transport  resistance  for  the 
substrate  due  to  turbulent  flow  and  diffusion  within  the  biofilm  and  biodegradation. 
The  turbulent  resistance  and  the  diffusion  term  are  both  approximated  by  Fick’s 
first  and  second  law  respectively.  Biodegradation  is  expressed  by  a  Monod  term, 
relating  substrate  consumption  to  a  reaction  rate  constant,  bacterial  mass  and  the 
half-velocity  constant  of  the  biochemical  reaction; 

dC  kXC 
at  ~  K.  +  C 
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Where: 


C  is  the  substrate  concentration,  \mgll\ 
fc  is  a  bacteria]  rate  constant,  [l/s] 

X  is  the  bacterial  density,  (mg//| 

Ke  is  the  half  velocity  constant,  \mgH] 

The  model  is  based  on  the  simplified  assumption  that  bacteria  attached  to  a  surface 
form  a  biofilm  of  a  homogeneous  thickness.  It  is  also  assumed  that  there  exists  a 
minimum  substrate  concentration  Cmin  below  which  no  biofiim  can  be  supported. 
A  mass  balance  considering  bacterial  growth  and  decay  allows  one  to  calculate  the 
biomass  from  the  substrate  concentration. 

The  effects  of  different  pore  sizes  of  porous  media,  different  flow  velocities,  and 
different  input  concentrations  were  simulated  with  the  model.  The  simulations  as¬ 
sume  a  steady  state  situation  of  substrate  concentration  where  the  transport  of  the 
substrate  is  in  equilibrium  with  its  bacterial  consumption.  The  results  are  shown  in 
Figures  15  to  17.  A  different  substrate  input  concentration  affects  the  substrate  and 
bacterial  distribution  within  the  column  significantly.  Variations  in  fluid  velocity 
and  grain  size  diameter,  however,  seem  to  have  a  less  pronounced  effect  on  the  sub¬ 
strate  concentration  distribution.  Depending  on  the  minimum  concentration  Cmin 
below  which  bacterial  growth  cannot  be  supported,  the  bacterial  mass  distribution 
differs  greatly  from  the  substrate  concentration  distribution.  Further  application  of 
the  model  is  planned  to  design  future  experiments  with  different  input  parameters. 

9.  Transfer  of  Laboratory  Study  to  a  Larger  Scale 

The  experimental  results  so  far  evaluated  apply  for  well-defined  laboratory 
conditions.  On  a  larger  scale  the  heterogeneities  and  complexities  of  an  aquifer  have 
to  be  taken  into  account.  Three  factors  that  are  important  here  include  the  average 
hydraulic  conductivity,  well  clogging,  and  multiple  electron  acceptor  conditions 
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Figure  15:  Biofilm  simulations  for  different  input  concentrations.  Substrate  con¬ 
centrations  and  biofilm  thicknesses  ewe  calculated  with  a  biofilm  model  for 
different  substrate  input  concentrations. 
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Figure  16:  Biofilm  simulations  for  different  specific  discharges.  Substrate  con¬ 
centrations  and  biofilm  thicknesses  are  calculated  with  a  biofilm  model  for 
different  specific  discharges. 
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Figure  17:  Biofilm  simulations  for  different  grain  diameters.  Substrate  concentra¬ 
tions  and  biofiim  thicknesses  are  calculated  with  a  biofilm  model  for  different 
grain  diameters. 
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Average  Hydraulic  Conductivity 

It  can  be  demonstrated  that  a  hydraulic  conductivity  reduction  on  a  small  scale 

•  of  the  order  of  millimeters  can  have  a  major  impact  on  a  larger  scale  of  the  order 
of  meters.  The  average  hydraulic  conductivity  of  layers  of  different  hydraulic  con- 

*  ductivities  arranged  in  series  perpendicular  to  flow  can  be  expressed  by  following 
equation(Freeze  and  Cherry,  1979): 

= ,4  w 

^  K, 

Where:  Kav  is  the  averaged  hydraulic  conductivity,  [m/s] 

c.  is  the  total  thickness  of  the  column,  [m] 

Ki  is  the  hydraulic  conductivity  for  layer  i,  [m/s] 

d,  is  the  thickness  of  layer  i,  [m| 

In  our  laboratory  experiment  with  a  column  of  10  cm  length,  the  hydraulic 
conductivity  was  reduced  by  six  orders  of  magnitude  in  the  first  cm  above  the  in¬ 
flow,  and  the  hydraulic  conductivity  of  the  other  9.0  cm  of  the  column  remained 
unchanged.  Applying  equation  (5)  for  the  average  hydraulic  conductivity  for  two 
layers  of  1  cm  and  9.0  cm  length  respectively,  one  can  calculate  a  reduction  of  the 
overall  hydraulic  conductivity  of  five  orders  of  magnitude  for  the  column  of  10  cm 
length  (Figure  18).  The  calculated  amd  the  measured  values  of  the  average  hy- 
,  draulic  conductivity  were  compared  for  the  experiment  with  an  input  concentration 

of  20  mg/1  acetate.  A  hydraulic  conductivity  decrease  by  six  orders  of  magnitude 
•  was  measured  in  the  first  column  segment  of  1  cm  length  while  for  the  total  length 

of  the  column  a  decrease  of  about  five  orders  of  magnitude  was  observed  (Figure 
19).  The  measured  hydraulic  conductivity  value  averaged  over  the  column  length 
corresponds  to  the  a  calculated  averaged  value  for  two  layers  consisting  of  the  first 
column  segment  and  the  rest  of  the  column.  If  the  same  reasoning  is  applied  for 
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Av«rage  Hydraulic  Conductivities 
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Figurt;  18:  Aver&ge  hydraulic  conductivity.  The  average  hydraulic  conductivity  of 
a  large  scale  can  be  largely  infiuenced  by  small  scale  hydraulic  conductivities. 
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Figure  19:  Hydraulic  conductivities  measured  over  time  for  the  first  column  seg¬ 
ment  (a)  and  for  the  toUl  column  length  (b)  An  input  concentration  20  mg/l 
acetate  is  represented  here. 
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columns  of  1  m,  10  m,  and  100  m  length  a  hydraulic  conductivity  reduction  of  six 
orders  of  magnitude  in  the  first  cm  of  the  column  induces  overall  hydraulic  conduc¬ 
tivity  reductions  of  four,  three,  and  two  orders  of  magnitude  respectively.  In  other 
words,  a  reduction  of  hydraulic  conductivity  at  a  small  scale  of  cm’s  can  have  a 
major  impact  on  the  hydraulic  conductivity  of  a  large  scale  of  hundreds  of  meters. 

Well  Clogging 

A  typical  biorestoration  scenario  will  consist  of  several  wells  injecting  nutrients 
to  enhance  bacterial  growth  (Figure  3).  Well  clogging  is  a  well  known  phenomena 
(van  Beck  and  Kooper,  1980;  Oberdorfer  and  Peterson,  1985;  Smith  and  Tuovinen, 
1985).  Since  laboratory  experiments  showed  that  permeability  reductions  occur 
very  close  to  the  nutrient  source  it  seems  likely  that  injection  wells  could  plug  and 
produce  only  a  small  zone  of  reduced  hydraulic  conductivities  (Figure  20a).  In  this 
case  the  hydraulic  conductivity  reductions  would  only  have  a  minor  effect  on  the 
fluid  and  contaminant  transport  in  the  aquifer.  To  prevent  clogging  it  is  common 
practice  to  add  oxygen  in  form  of  toxic  hydrogen  peroxide  (Yaniga  and  Smith,  1984) 
or  ozone  (Nagel  et  al.,  1982)  which  prevents  bacterial  growth  close  to  the  well. 
In  travelling  through  the  aquifer  these  toxic  compounds  are  degraded  amd  oxygen 
is  finally  available  in  its  non-toxic  form  for  bacteria.  A  zone  of  lower  hydraulic 
conductivity  around  the  wells  might  form  Md  deviate  the  groundwater  flow  around 
the  well  (Figure  20b).  It  is  however  very  diflicult  *o  estimate  the  location  and  the 
amount  of  hydraulic  conductivity  reduction.  The  highly  oxidizing  compounds  will 
react  with  the  aquifer  matrix.  Hydraulic  conductivity  reduction  due  to  mineral 
dissolution  might  overlap  with  the  biological  effects.  The  transport  distance  of 
the  toxic  compound  will  depend  on  the  mineral  composition  of  the  aquifer,  the 
reaction  velocity  of  the  compound,  the  advective  and  dispersive  transport,  and  the 


Figure  20:  Possible  flow  patterns  around  injection  wells,  (a)  Permeability  reduc¬ 
tions  very  close  to  the  well  does  not  influence  the  Sow  pattern,  (b)  Due  to  the 
toxicity  of  hydrogen  peroxide  bacterial  growth  can  develop  only  in  a  distance 
r  from  the  injection  well.  In  an  extreme  case,  high  permeability  reduction 
could  deviate  the  polluted  groundwater,  preventing  its  degradation. 


adsorption  coefficient.  Bacterial  densities  are  difficult  to  define  in  the  transition  zone 
between  toxic  and  non-toxic  area.  The  biologically  active  zone  might  be  smeared 
out  due  to  diffusion  and  toxic  activity  and  may  therefore  have  a  smaller  effect  on 
hydraulic  conductivity,  if  at  *11. 

Multiple  Electron  Acceptor  Conditions 

The  above  described  laboratory  experiments  were  performed  under  the  assump¬ 
tions  that  aerobic  conditions,  i.e.  oxygen  served  as  the  electron  acceptor,  are  pre- 
vailent  in  a  large  scale  aquifer.  On  a  large  scale  the  situation  might  be  more  complex 
however.  Depending  on  the  background  concentration  and  the  supply  of  oxygen  as 
well  as  the  degree  of  bacterial  activity  in  an  aquifer,  oxygen  might  become  de¬ 
pleted  and  other  compounds,  such  as  nitrate,  sulfate  and  carbon  dioxide  take  over 
as  electron  acceptor  (Figure  7).  Thermodynamic  calculations  (Appendix  1)  as  well 
as  observations  from  field  experiments  (Oberdorfer  and  Peterson,  1985)  show  that 
under  denitrifying  conditions  (nitrate  as  electron  acceptor)  nitrogen  gas  and  un¬ 
der  methanogeneris  (carbon  dioxide  as  electron  acceptor)  methane  is  likely  to  be 
formed.  Permeability  is  reported  to  be  considerably  reduced  by  gases  (Christiansen, 
1944;  Gupta  and  Swartzendruber,  1962),  so  that,  both,  biomass  and  gases  might 
contribute  to  significant  permeability  reductions.  However,  it  is  hardly  possible  to 
predict  the  amount  and  extent  of  gases  produced  as  well  as  it  is  difficult  to  assess 
its  relative  importance  for  clogging  effects  compared  to  bacterial  mass  formation. 

Summary 

As  pointed  out  at  the  beginning  of  this  chapter  small  scale  permeability  re¬ 
ductions  can  have  a  major  impact  on  large  scale  permeabilties.  The  situation  in 
a  large  scale  aquifer  system  is  however  so  complex  that  it  seems  inappropriate  to 
quantify  effects  of  bacterial  clogging  in  large  scale  aquifers  at  this  stage  of  research. 
Despite  this  fact,  it  appears  neccessawy  for  any  design  of  biorestoration  methods 
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to  take  into  consideration  permeability  reductions  due  to  bacterial  activity  for  the 
following  reasons: 


•  It  is  very  likely  that  major  permeability  reductions  due  to  bacteria  actually 
occur  also  on  a  large  scale. 

•  Permeability  reductions  on  a  large  scale  would  alter  fluid  flow  and  therefore 
contaminant  transport  significantly.  Ignoring  permeability  reductions  could 
lead  to  erroneous  predictions  of  fluid  and  contaminant  transport. 

•  Although  quantitative  predictions  of  flow  and  contaminant  transport  are  not 
yet  possible  for  aquifers  under  the  above  described  conditions,  a  qualitative 
analysis  of  the  impact  of  permeability  reductions  due  to  bacteria  should  be 
mandatory  for  biorestoration  design. 

•  A  numerical  model  simulating  fluid  and  contaminant  transport  considering 
biological  clogging  ,  would  represent  a  useful  tool  for  biorestoration  design. 

A  first  attempt  at  numerical  modeling  is  described  below. 

10.  Large  Scale  Modeling 

Various  numerical  models  exist  that  include  biodegradation  in  the  contami¬ 
nant  transport  algorithm  (Sykes,  et  al.,  1982,  Borden  et  al.,  1984,  Corapcioglu  and 
Haridas,  1985,  Benefield  and  Molz,  1985,  Molz,  1986,  Borden  and  Bedient,  1986). 
No  attempt  has  been  made  so  far  to  include  hydraulic  conductivity  reductions  due 
to  bacterial  growth  in  any  numerical  model.  A  numerical  model  was  developed  here 
to  simulate  fluid  and  contaminant  transport  including  biodegradation  and  chang¬ 
ing  hydraulic  conductivities  due  to  bacterial  activity  (Appendix  1).  The  following 
differential  equations  describe  the  physical  and  biochemical  system: 
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The  balance  of  mass  equation  for  fluid  flow  in  two  dimensions  at  steady  state 
considering  changing  hydraulic  conductivities; 

where:  K{x,y,i)  is  the  hydraulic  conductivity,(m/sl 

h  is  the  hydraulic  head,  [m] 

The  balance  of  mass  for  the  contaminant : 


dc 

dt 


3  dC  „  dC,  d  dC  „  dC^  dC 
~  dx  dy  ^  dx^^'^dy^~  dx  ~ 


(7) 


dy 


where:  C  is  the  contaminant  concentration  ,  [my//] 

Dii  is  the  dispersion  coefficient  tensor,  [m*/s] 

k  is  the  rate  constant  for  biodegration,  [l/s] 

u,  and  V,  is  the  average  linear  velocity  components,  [m/s] 


K(x,y^t]  dh 
n  dx 


(8) 


(9) 

n  dy  '  ' 

where:  n  is  the  porosity  of  the  porous  medium,  [-] 

The  steady  state  form  of  (6)  is  based  on  the  assumption  that  the  hydraulic  head 
equilibrates  very  quickly  to  changes  in  permeability,  on  a  tune  scale  much  shorter 
than  the  solute  transport  in  equation  (7).  The  linearity  of  the  biodegradation  term 
is  justified  if  dealing  with  relatively  low  concentrations.  The  Monod  expression 
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(10) 


dC  k!XC 
dt  ~  Kc  +  C 

where:  k'  is  the  max.spezific  rate  of  substrate  utilization,  [l/s] 

X  is  the  bzLcterial  density,  (mg//) 

Kc  is  the  Monod  half  velocity  constant,  (my//) 
normally  used  to  describe  biodegradation,  becomes  approximately  linear  if  C  >  Kc, 
i.e.: 


dt  ~  Kc^ 


(11) 


where: 


(12) 


Depending  on  the  level  of  biological  activity  the  hydraulic  conductivities  can  be 
either  constant  or  change  over  time.  The  following  equation  approximately  describes 
the  exponential  decrease  of  hydraulic  conductivity  being  observed  in  our  laboratory 
experiments; 

K  =  (13) 

Where:  Kq  is  the  initial  hydraulic  conductivity,  [m/s] 

A  is  a  rate  constant  for  hydraulic  conductivity  reduction,  [l/sj 

A  Galerkin  finite  element  model  was  developed  to  solve  equations  (6),  (7),  (8), 
(9),  and  (13).  This  model  was  applied  to  test  the  significance  of  hydraulic  conductiv¬ 
ity  reductions.  A  fictitious  aquifer  was  modeled  Including  a  contamination  source 
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and  an  area  of  increased  biological  activity  with  reduced  hydraulic  conductivity 
values  (Figure  21). 

The  contaminant  leaking  from  its  source  was  transported  with  the  groundwater 
and  eventually  reached  the  biologically  active  zone.  It  was  assumed  that  this  zone 
was  activated  through  infiltration  of  nutrients  and  that  the  infiltration  rate  was  so 
low  that  it  did  not  affect  the  fiuid  fiow  system  in  the  aquifer.  With  infiltration  the 
hydraulic  conductivity  decreased  exponentially  with  time.  One  example  was  calcu¬ 
lated  where  infiltration  of  nutrients  was  initiated  150  days  after  the  contamination 
occurred.  Biodegradation  was  not  explicitly  considered  in  this  example.  Com- 
parision  of  the  contaminant  transport  of  both  cases,  with  and  without  biological 
activity,  show  that  the  flow  field  is  significantly  altered  (Figure  22)  and  the  contam¬ 
inant  plume  is  retarded  and  dispersed  after  500  days  of  contamination  (Figure  23 
and  Figure  24). 

Different  scenarios  and  their  effect  on  contaminant  transport  would  have  to  be 
simulated  in  future  studies  because  of  the  uncertainties  in  quantifying  the  hydraulic 
conductivity  reductions  without  field  data. 

11.  Conclusions 

Bioreclamation  of  polluted  groundwater  ought  to  be  a  successful!  and  cost- 
effective  treatment  method.  An  overview  of  common  technologies  for  in-situ  treat¬ 
ment  demonstrates  that  the  success  of  the  method  will  be  closely  linked  to  the 
hydrogeological  conditions  of  the  polluted  aquifer.  One  purpose  of  this  work  was 
to  demonstrate  that  bioreclamation  schemes  czm  be  significantly  affected  through 
permeability  reductions  induced  by  bacteria.  Biochemical  reactions  in  the  subsur¬ 
face  induced  through  the  infiltration  of  nutrients  for  the  enhancement  of  bacterial 
growth  can  alter  hydrogeologic  parameters,  especially  permeability.  A  literature 
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Figure  21:  Finite  element  mesh  for  model  aquifer  with  contairunant  source  (red) 


and  biodegradation  area  (green). 
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Figure  22;  Numerical  model  results:  (a)  hydraulic  h^ad  pattern  and  (b)  velocity 
field  500  days  after  contamination  and  350  days  after  biological  growth. 
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Figure  23:  Spread  of  contamination  250  days  after  pollution  incident:  (a)  with 
bacterial  growth,  (b)  without  bacteria!  growth 
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Figure  24:  Spread  of  contamination  500  days  after  pollution  incident:  (a)  with 
bacterial  growth,  (b)  without  bacterial  growth. 
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review  has  demonstrated  that  permeabilities  can  be  reduced  by  several  orders  of 
magnitude  due  to  bacterial  growth.  In  laboratory  experiments  we  confirmed  this 
observation.  Most  of  the  permeability  changes  occurred  very  close  to  the  nutri¬ 
ent  source.  After  a  lag-phase,  the  permeability  appeared  to  decrease  exponentially 
with  time.  The  substrate  input  concentration  has  a  large  influence  on  the  degree 
and  rate  of  permeability  reductions.  Further  experiments  are  needed  to  confirm 
the  observations  made  so  far.  The  influence  of  other  par2uneters,  such  as  grain-size 
and  fluid  velocity  on  permeability  changes  need  to  be  identified.  The  influence  of 
fractures  is  unknown. 

Difficulties  are  likely  to  be  encountered  in  transfering  the  small  scale  results  into 
the  larger  scale  of  an  aquifer.  The  geological  and  biological  systems  in  aquifers  are 
generally  so  complex  that  it  seems  impossible  at  this  stage  of  research  to  apply 
the  laboratory  data  in  the  field  without  further  assumptions.  Laboratory  and  field 
tests,  however,  suggest  that  permeability  reductions  actually  occur  in  the  subsur¬ 
face.  Although  quantification  at  this  moment  is  not  possible,  a  consideration  of  the 
effects  of  permeability  reductions  due  to  bacteria  seems  to  be  neccessary  for  any 
description  of  contaminant  transport  during  biorestoration  scenarios.  A  numerical 
model  with  the  ability  to  test  different  assumptions  concerning  amount  and  ex¬ 
tent  of  permeability  reductions  would  represent  a  useful  tool  for  the  design  of  any 
biorestoration  scheme.  Field  experiments  could  be  used  to  validate  the  model  for 
predictive  purposes. 

The  here  described  research  is  to  be  considered  as  one  step  towards  a  better 
understanding  for  the  overall  hydrologic  system.  Further  laboratory  and  field  ex¬ 
periments  are  needed  in  order  to  fully  understand  the  role  of  bacterial  clogging 
effects  on  a  large  scale.  Regardless  of  the  scientific  value  of  quantifying  this  phe¬ 
nomenon,  the  hydrogeologic  effect  of  permeability  reduction  will  undoubtedly  play  a 
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significant  role  in  the  practical  application  of  in-situ  biorestoration  of  contaminated 
groundwater. 
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Transport  model  program  listing  and  sample  input  data  deck 


62 


M  M 
M  M  M 


VI  M  M 
•I  M  VI 


a 


o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  e 
o  o  o 


e  o  o 
o  o  o 
o  o  o 
o  o  o 


o  o  o 
e  o  o 
o  o  o 
o  o  o 
o  o 


o  •  o 
o  o 
o  >  o 
o  o 
o  fi  e 
o  o  e 


0  9  0 
O  >  O 

o  o 
O  VI  o 
e  X  e 
o  >  o 
0X0 
o  M  o 
o  «<  o 
o  >  o 


»• 


M  M 

M  M  M 

U  U  (J 

MUM 
(M  M  M  U  M  M  U 


o  a  o 
o  o  o 

O  '<4  O 


o  o  o 

O  «  0  A  O  e  A 


Cl  u 

VI  u 

U  Cl 

Cl  Cl 

u  u 

U  Cl 

VI  u 

U  Cl 

VI  Cl 

U  Cl 

U  Cl 

u  u 

u  u 

u  u 

U  VI 

VI  u 

uuuuuuouou 

uuuouuuuuu 


#1  ^  0  0 


>90099999999 

>99999999999 


9  9 
9  9 
9  9 
9  9 
9  9 
9  9 
9  9 
9  9 
9  9 
9  9 


0009000999 

0000900000 


O  O  O  fl 

O  O  O  C 

o  o  -  o  c 

o  O  0  9 

9  0  O  O 

oooooooooooooo 

oooooooooooooo 


A  A 
A  A 
A  A 
A  A 
A  A 
A  A 
A  A 
A  A 


AAA 

AAA 

A  A 
A  A 
A  A 
A  A 
A  A 
A  A 


A  A 
A  A 


A 

A 


A 

A 


A 

A  A 
A  A 


A  V  A 
AAA 


o  o  o 
O  VI  o 


0  6  0 
O  •  o 
o  A  O 
O  OO 

O  o 
0  9  0 
o  c^o 

O  M  O 
O  O 
0-4  0 
O  A  O 
O  44  O 
O  *4  o 
o  oo 
O  •i*  o 
o  O  o 
o  o 
o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 
o  o  o 


<4  M 

I  M  M  U  M 


AABASAAXBBBBZZ 

SSXASABAZZABZK 
X  X 
K  X 

X  X 

m  m 

X  K 
X  K 

XXXXXXXXXKXXXX 

XXAKZXXXKXXXKX 


o  o 
o  o 
o  o 
o  o 
o  o 
o  o 


oooooooo 

oooooooo 


oooooooo 

oooooooo 


o  o 
o  o 


o 

o 

o 

o 

o 

o 


o  o 
o  o 


o 

o 

o 

o 

o 

o 


VI  VI  VI 

A  0  A 

AAA 
AAA 
A  A  A 


»  X  »  V  »  S  » 


o  o 
o  o 
o  o 
o  o 
o  o 
o  o 


oooooooooo 

oooooooooo 


oooooooooo 

oooooooooo 


o  o 
o  o 
o  o 
o  o 
o  o 
o  o 


o  u 
u  u 
u  u 
u  u 
u  u 

VI  VI 
VI  VI 
U  VI 


vr 

v; 


VI  VI 
VI  VF 


VI  VI  VI  VI  VI  VI  u 
U  U  VI  VI  U  U  VI 


V  u 
VI  u 
VI  Cl 
VI  VI 
VI  u 

V  o 
VI  u 
U  V 


A  k. 


AAA 

AAA 

AAA 

AAA 

AAA 

AAA 

AAA 

AAA 

AAA 

AAA 


AAA 

AAA 

AAA 

u 

AAA 

M 

AAA 

9 

>1 

AAA 

A 

AAA 

K 

A 

AAA 

Jt 

to 

AAA 

• 

AAA 

A 

•  • 

to 

0  m 

Oi 

»  • 

A 

e  O  O 

0 

o  e  o 

•  J» 

e  e  o 

04 

B 

e  o  o 

•A  A 

u 

e  o  o 

♦*  0 

B 

o  o  o 

B 

O  A  o 

04  m 

M 

o  e  o 

m  •A 

> 

O  O  0 

•A 

o  o  o 

A 

o  o  o 

0  A 

S 

e  o  o 

•  o 

3 

o  o  o 

9  U 

O 

e  o  o 

9  • 

U 

e  o  o 

•  A 

u 

o  o  o 

• 

0 

o  o  o 

o  e  o 

M  m 

to 

o  o  o 

u  • 

0 

e  o  o 

o  «• 

o  o  e 

•4  a 

0 

o  o  o 

o  o 

9 

o  o  o 

04 

o  o  o 

0 

b 

o  e 

p*  • 

— 

e  o 

»4 

O  0  o 

•  H 

B 

o  •  o 

M 

O  40  o 

• 

B 

O  >  0 

•A  • 

U 

0  O 

04 

A 

O  fi  O 

•  • 

A 

0  0  O 

•A  W 

» 

0  1*4  0 

0  A 

* 

0  0  0 

••  s 

B 

0  to  o 

o  e 

U 

0  0  O 

04  U: 

> 

0  >  0 

B 

0  O 

• 

< 

O  A  O 

A  <9^ 

O 

O  X  O 

A  « 

o 

0  >  O 

A  'A 

0  \  O 

1  A 

O  M  O 

X  U 

u 

o  4  e 

9  • 

A 

o  >  o 

n  o 

1 

9 

o  o 

O  1  o 

A 

• 

O  0 

A  « 

H 

0  0  0 

u 

0  • 

0  O  O 

6 

M  0 

O  ■«4  O 

0 

0  < 

O  44  O 

M  M 

0  0  0 

^  • 

A  f* 

e  to  o 

•  A 

A 

0  0  o 

•  A 

>  • 

O  BO 

A  Ck 

9 

O  to  o 

>  ■ 

to  • 

O  0  O 

•  A 

•  9 

O  VI  o 

•  9 

o  o 

IS 

9 

0  41  o 

A  A 

■ 

0  6  0 

•  A 

>«  0 

O  0  O 

•  > 

J>  to 

O  fi  o 

04 

M 

O  BO 

M 

A 

O  •«4  O 

»  A 

0  A 

0  9  0 

0* 

•  •  0 

O  90 

e  c 

e  *• 

O  U  o 

•  • 

A  O 

0  O 

0  A 

A 

0  o 

• 

r* 

O  0  o 

0  r* 

O  44  O 

A  A' 

0  0 

O  >4  o 

A  9 

A  0 

O  90 

0  9 

1  A 

O  *4  o 

»■»  A- 

B  1 

o  o  o 

te 

9  B 

0  O 

0  • 

1^  9 

o  o 

••  > 

1  0 

O  O  0 

m 

0  1 

o  o  o 

O  9 

A  0 

o  o  o 

lb  W 

A 

o  o  o 

•  • 

0 

o  o  o 

u 

O  6 

o  o  o 

9  A 

0 

o  o  o 

O  7 

0 

o  o  o 

B  %0 

4  •• 

o  o  o 

O  Q 

M  0 

o  o  o 

Cl  U 

»•  4 

o  o  o 

—  • 

M 

o  o  o 

5  *• 

O  |4 

o  o  o 

4 

f 

O  0  o 

«  • 

o  o  o 

^  k4 

o  o  o 

O 

•  0 

o  o  o 

m 

9  A 

o  e  o 

•  0 

0  o  o 

•  • 

9  A 

o  o 

m  — 

»* 

o  o  o 

M 

o  o  o 

B  m 

o  o  o 

m  u 

A  e 

o  o  o 

U  0 

0  9 

>  M 

A 

•—  M 

•« 

•'  5 

• 

^  9 

U  4i 

AAA 

4  • 

9  u 

AAA 

9  B 

Q  0 

AAA 

o  w 

B  A 

AAA 

i> 

O  0 

AAA 

m 

U 

AAA 

9  4 

• 

AAA 

A  o 

O  o 

AAA 

-4  O 

o  o 

AAA 

B  — 

1^  A 

AAA 

UC  uses  THE  riaXTC  SLeMEMT  METHOD.  TRXA.HaULXR  CLEHCRTS 
LINEAR  INTERPOLATION. 


I 


f 


I 

I 


uuuuuouuuoououuoo 


I 


/ 


u 

< 

H 

V) 

>• 

o 

< 

u 


o 


0X^30 

OOF*  o  m  Q 

at  X  M  MO  OB 

M  4  N  At  (/)  U  4 


u 


« 

« 

X 

M 

>• 

« 

M 

O 

4 

H 

«  «  • 

F* 

o 

X  . 

• 

* 

w 

«  .»  A 

o 

>•  ♦ 

H 

V) 

«  *«  o 

% 

O 

VI 

B 

m 

«  « o 

O 

•  o 

O 

O 

M 

«  o  o 

O 

8 

M  4 

M 

M 

> 

«  o  ^ 

8 

9 

4  H 

O 

M 

•  O'— 

8 

X  4 

8 

4 

8 

«  OK 

» 

8 

4 

9 

9 

«  w  t« 

M 

« 

o 

k 

o» 

«  u  u 

O 

8 

%  ^ 

8 

4 

u 

V) 

«  •  8 

o 

W  O 

4 

M 

8 

•  —  > 

w 

« 

H  O 

B 

O 

H 

M 

«  <« 

8 

M  O 

4 

8 

4 

«  «  • 

O 

o 

8  ^ 

VI 

U 

Oi 

•  O  — 

u 

o 

>-  — 

A 

k 

X 

O 

«  o  o 

X 

a 

O 

— > 

V  •  » 

o 

k 

4 

u 

8 

«  o  o 

* 

M 

•  4 

o 

4  8  4 

8 

o 

O 

•J 

«  r4  o 

H 

M  F< 

o 

4  9  9 

4 

— * 

8 

B 

w 

(/) 

1^ 

8 

M  4 

o 

0  4  0 

M 

iJ 

O 

X 

8 

•  H  — 

« 

O 

M 

•  •  « 

X 

U} 

o 

4 

w 

X 

•  X  >•  N 

8 

M 

— • 

8  • 

— 

OK  OK  01 

9 

4 

H 

O 

•  •  9  M 

M 

O  o 

4  — 

4 

8  8  8 

» 

8 

k 

M 

«  A  • 

8 

o 

o 

u 

4  4  4 

4 

»A 

k 

B 

*  O  -» 

* 

8 

o 

«  VI 

01 

»  M  W 

8 

O 

» 

M 

4 

M 

• 

•  O  O 

8 

« 

« 

4 

N  ■  a 

4 

O 

U 

> 

H 

•m 

H 

«  O  O 

o 

w 

8 

4  •  « 

8 

M 

8 

8 

8 

M 

«  r<  O 

o 

* 

O 

F* 

M  *1 

H 

• 

WWW 

O 

M 

U 

e 

O 

o 

4 

8  8 

8 

— • 

rw  rf» 

« 

% 

U 

k 

k 

X 

M 

•  X  — ’ 

u 

o 

X 

O  *-• 

M 

o 

0  O  0 

•J 

> 

8 

U 

4 

t« 

W 

«  *4 

w 

X 

8 

o 

8 

o 

WWW 

U 

8 

a 

X 

H 

«  —  o»  U 

M 

U 

« 

O 

o 

a  8  8 

> 

« 

o 

8 

o 

VI 

«  o  • 

M 

F* 

m 

4 

•  4 

u 

•  •  • 

8 

« 

8 

u 

*4 

H 

•  o  — 

a 

8 

o 

M  «J 

WWW 

X 

u 

8 

4 

m 

«  e  o 

o 

H 

• 

8  U 

4 

•  •  • 

» 

8 

u 

O 

M 

«  o 

M 

• 

8 

O  O 

» 

U 

e  fl  a 

V 

a 

k 

k 

M 

8 

•  —  o 

o 

U 

4 

H 

WWW 

8 

« 

o 

8 

kJ 

V) 

M 

•  8  V 

• 

y 

> 

o 

• 

8 

• 

U 

•  •  • 

O 

X 

u 

4 

X 

9 

M 

•  •  >» 

Q 

•• 

o 

«  4 

O 

— • 

M 

>  »  > 

{J 

u 

•J 

M 

8 

N 

•  —  ►*  O 

O 

8  iJ 

U 

e 

8 

w  ^  mm 

8 

M 

k 

4 

kJ 

VI 

«  o  > 

o 

• 

o 

O  4 

u 

o 

1 

M 

> 

a 

8 

4 

M 

H 

*  o  • 

o 

o 

U 

o 

U  O 

o 

F* 

^ 

8 

« 

8 

8 

8 

«  p4 

o 

4 

4  M 

« 

9 

«  •  4 

o 

8 

« 

k 

k 

8 

*  *0 

o 

H 

a 

•w 

X 

8 

8 

9  9  9 

M 

4 

k 

» 

M 

X 

9 

O 

8 

«  o  o 

8 

n 

M 

8 

8 

8 

•  8 

O 

fal 

F» 

^  ^  ^ 

4 

« 

U 

U 

O 

M 

8 

H 

U 

«  o  o 

u 

X 

o 

O 

hI 

U  M 

V 

a 

9 

*  «  » 

8 

4 

o 

8 

o 

U 

8 

H 

H 

•  o  ^ 

u 

F* 

H 

u 

u 

4 

8  X  U 

\ 

O 

%  «  « 

o 

8 

o 

O 

8 

U 

H 

4 

H 

«  r«  «- 

8 

J  S»- 

N 

H  4 

M  4 

41  *4  r* 

w 

8 

U 

a 

8 

o 

8 

M 

•  8 

O 

4  w 

X 

N 

o 

o 

^  -w  w 

W 

8 

8 

M 

8 

«  4  > 

o 

X 

4 

»4 

4 

8 

8  8  8 

8 

8 

•«> 

•V 

X 

> 

M 

« 

O 

w 

4 

O  O 

M 

O 

a 

X 

M 

m 

«  8 

F* 

4 

iJ 

H 

M  M  M 

4 

M 

O 

4 

8 

4 

«J 

VI  9t 

«  O 

U 

8 

8 

8 

1 

8  8  8 

9 

4 

M 

O 

9 

« 

ii 

4  •H 

«  M 

s 

N. 

s 

N 

»« 

8  8  8 

O 

8 

01 

U 

O 

4 

» 

«  VI 

8 

8 

8 

8 

9 

4  4  4 

OK 

O 

(9 

8 

a 

8 

B 

«  8 

O 

O 

O 

O 

8 

m 

4 

X  'O 

«  M 

X 

X 

X 

X 

8 

k4  kJ  U 

9 

•J 

U 

•J 

iJ 

M 

«  X 

X 

X 

X 

X 

M 

U  «J  U 

8 

•J 

•J 

4 

4 

4 

V) 

« 

«  M 

O 

o 

o 

O 

4  4  4 

4 

4 

4 

4 

4 

«) 

O  M 

«  o 

u 

u 

u 

u 

.4 

u  u  o 

ij 

O 

U 

U 

u 

O 

9 

O  8 

«  • 

• 

« 

m 

• 

• 

« 

o 

• 

4 

M  9 

4 

4 

O 

Oi  ») 

U 

u 

ouuooouuuuouuuuuu 


O  V  V 


OOl.;  UUU  UUU  UUU  UOU  L»UU 


OPTION  OP  READING  IN  VELOCITIES 


CALL  VEL  (HELEK*  VX ,  VY .  V8ETB) 
CO  TO  I 


V  U  U 


o  u  u 


{j  u  u  uuo  oou  ouu  you  uuu  uoo  uuu 


u  u  u 


uuu 


It  COHTINUE 
14  COSTXHUE 


COHTIMUt 


OPTION  IP  SOLUTE  TEAMSPOBT  OB  NOT:  laPLUID  TBANSPOBT  OHLT 

2-SOLUTC  TBANSPOBT 


IS  All  CXrOHIHTIAL  MtYbSAULIC 
DECkEASC  AT  TIMISTSP  RTH 


cooBomATEs  rot  becuus  cbid 


TO  SOS 


SSS  COflTXHUC 


>EAD  BLCHSMT  IVCIDMCES 


COHCK1IT&ATZOH 


I 


ICOOCH(M)>KOOlHt 
H(N)-UEAO 
11  COHTXBUI 


»• 

m 

o 

a, 

VI 

z 


H 

9 

O 


> 

o 

z 


>• 

M 


H 

9 


m 

< 

Q 

Z 

9 

O 


9 

tJ 


>4 

9 


U 

o 

o 


M  #4 
9  • 


l«  9 
Z  < 
O  M 

U  « 


^  U  *4  w  r4  #»• 

•OHX<50ai3 
'^Z99ZZZZ 
^  IkikJWMUM 

355§;§i§s 

Q*4f«*,4^f^»r4Z 

OOUUUUUUO 

OMOOMHMNU 


o 

t« 


o 

9 


»•  »«  B 

••  «  Z  H  O 

jS  P  p  M  B 

m 

o 

lb 


z  k.  B  I 


U  O  U  V  U  U 


*  o 

>  m 

•  o 

'  z 

4 

m 

e 

z 

w 


H 

S  ^ 


z 

M 

H 


9 
*  9 


2  2  “ 


I  .40 
>  O  O 
i  B  tt  • 
'  O  Z  ^ 
U  4  o 
Z  O 
•  o 
J  •  ^ 

M  M 
>  Z  V) 

■  z  o 


w  ’ 
Z  ' 


9 

*  4 
I  H 
t  W 


«  u  z  z  e 


e  M  • 


I  o  z  • 
>  o  u 

•  B  4  O 

o 

*  •  o 
B  J  V 
4  9'-' 
«  z  « 

Z  9  o 

B  4  Z 


z 

B 


5  :: 


z  9  • 
M  Z  -> 
U  Z  O 
It  M  O 
BOO 


B 

4 


9 

or 

4 


9 

O 


9 

U 

s 

5 

B 

4 

Z 

N 

B 

O 

5 

o 


^  Z  p« 

w  B  w 

Z  M  O  I 
4  B  Z  I 
B  O  M  I 
Z  U  I 
J  •  ( 
4  «  Z  < 
M  O  I 
«  B  M 
—  90 
O  9 
o  * 
o  '  M 
^  B  iJ 
O  B 
^  O  O 
4  M 
B 

z  %B 
O  M  M 

4  Z  Z 

M  Z 

N 

Z 

O 

iJ 

z 

N 


O 

u 

z 


z 

«l 


:  V 

•  u 

I  z 


B  Z 
B  O 


o  o 


o 

z 


B  Z  O  O 
O  M  B  9 

U  S  O 
o  z  u  ^ 
B  NZ  -4 


V 

\ 


o 

s 

o 

u 


z 

z 


z 

9 

9 

O 

B 

O 

U 


kJ 

9 

4 

B 

Q 


O 

4 


9 

9 

4 

Z 

H 


O 

o 


4 

Q  n 


►I  f4 


z 

B 

H 


o  z 

z  o  z 

Z  M 

o  o 

9  0  U 

9 


9 
9 

2  ^'0  4 

O  o  w  B  Z 

9  O  ri  O  W 

m  M  H  I  O  I 

^9  -Z  ■  - 

B  tfi  B 

^  M  w  o  Z  ^ 

Q  Z  O  ^  ^  ^ 

4  B  4  B  4 

M  O  M  O  O  Z 

«  U  B  O  U  W 


z 

M  ^ 
O  • 

M  9 
B  M 

Z  B 


4 

B 

O 

Z 


f4 


U  U  U  U  U  9 


9  9  U 


COMTIMUE 


KEAO  ZM  rOROSITT,  DXPPUSIOM  CORPPZCIINT,  LONaZTUOZMXL  DISPERSZVZTY 
TRANSVERS  DZSPERSIVITY 


CREATES  A-MATRZX,  Q>VCCTOR,  CALLS  HATRZX  SOLVER,  TRAHSPCRS  RESULTS 


MOO  • 

o  o  o 

o  o  o  o 

o  ^  ^  o 

«  V  w  O 

t  ^  k  « 

«  B  y  u 

n  a  Q  VI 

-  o  o  o 

;  «  w  M  «  ^ 

J  M  O  o 

)  O  •  «  flh  o 

;  o  —  *-  o 

I  o  o  o  •  « 

^oo—  — 

*  O  O  O  H 

•  B  ^  ^  O  ^ 

n  M  ^  ^  O  S 

9  a  <«>  a. 

»  •  y  U  w  »j 

I  ^  Q  O  «  •< 

» o  o  o  o 

O  M  M  H  • 


W  M  O 

•N  o  o 

>•  a  ^ 

o  ^ 

at 

• - H 

•  '-x  o  X  y 
M  o  o  X 
o  o  o  X  > 

o  o  »  -X 

o  ^  o  - 

•-x  r<<  y  o  o 

9  y  X  <*«  o 

-<  X-  o  - 


- -  y  » 

O  O  O  >  M 
O  O  O  O 

O  O  O  O  •  X 

^  xp  ^  ^  M  M 

W  ^  .  X  O  O 

r«  .-4  y  o  o  o 

9  H  O  r»  a  o 

o  ^  ^ 


BXr«x*  O  0^9 

yH'-'x-x'^'O  .IN  —  w 

tJy^OOx-HO  •  xM  —  yy 

bia  ooo'«-xx-xo<<>9 
M>  xaoox-oooo 
—  -*XXPN*0<JOOO  X  •  X 

Ox-.— .^4oo«r  —  —  -. 

y  0XUx-.9xp«— ooo 

a  oyyBA>x--.Mooo 


J  a  hJ 

y  y  w 

y  M 

—  w  —  y 

M  M  M  B 


< 

IN  Q 

a 

O  iJ 

•H  ^ 

B  o 

o 

o 

M 

X 

•4  U 

.J  U 

xJ 

-X  H 

y 

^  O 

o 

u  -< 

9  M 

N» 

y 

o 

M 

M 

w  y 

y 

y  y 

y 

H 

M  W 

y 

\ 

\ 

y  X- 

X— 

— 

o 

8 

8 

•J 

n  8 

M 

8  M 

y 

y 

X  y 

M 

v.\ 

B 

a 

O  M 

M 

O 

4  O 

H 

M 

o 

V 

#4 

M  8 

y  M 

8 

B 

4  y 

ij 

9 

8 

8  ' 

4 

8 

y  y 

y 

M  M 

M 

—  — 

tf 

H  O 

8 

4 

y 

M 

8 

M  O 

i 

. 

4 

•  > 

4  • 

8 

y  o 

4 

8 

H  a 

M 

X 

. 

•4 

M 

M  M 

f  ■ 

■ 

■  • 

• 

y 

• 

•J  o 

iJ 

•>4 

M  H 

8  M 

X 

8  O 

8 

*4  1 

w 

M 

$ 

XXX 

8 

y  y 

to)  1 

M 

1 

y 

»« 

S.  % 

\ 

\ 

s  *-• 

»« 

• 

f 

9 

u 

X 

— > 

M  M 

M 

M  M 

M 

4 

M 

8 

M 

•J 

9 

8  8 

8 

B 

B  8 

4 

M 

PJ  M 

a 

8 

8 

8  8 

xJ  iJ 

x4 

«J  u 

o 

X 

1  ■ 

t 

8 

8  ^ 

M 

8 

O 

O  O 

O 

O 

O  8 

X 

►1 

M 

4 

O 

M 

1^  8 

M  y 

y 

y  y 

y 

y 

M 

8 

►I 

M 

8 

5  5 

X 

X 

X  M 

B 

o  - 

y 

-« 

1 

•  • 

M  n 

8 

►4  1^ 

8 

y 

4 

^  04 

m 

04  8 

04 

8 

M 

. 

y 

SC  X 

X 

X 

8  X 

4 

M 

8 

U 

iJ 

•J  U 

X—  <— 

w 

—  — > 

— 

4 

y 

— .  — 

— 

M 

• 

M 

8 

9 

o  o 

o 

o 

O  IN 

o 

O  'x- 

O 

9 

o  y 

M  y 

y  y 

y 

y  y 

y 

9 

8 

8  8  8 

o 

0  8  1 

y 

— > 

O 

VI 

u  u 

fj 

u  o 

O 

0  0  4 

U 

O 

o 

M 

n  8 

m  m 

8 

u  u 

u 

•J 

4 

M  M 

M 

o 

M  O 

M  n 

X 

4 

u 

• 

0 

« 

8 

8 

4 

« 

M 

m 

O 

> 

o 

M 

O 

y 

to4 

COMDUCTAHCK  HATBIX  AT  COHSTABT  HEAD  BOUNOABT  COHDITIOMS 
(DiaiCHLET  BOUBOABT  COHDITtOHS) 


35  COMTIKUC 


SPCCiriC  DISCKAKSE 


ABIA(N)  -  (CC(  JED  •■KICL)  -  CE  ( I  EL)  *EE  (  JED  1/2 


u 

o 

•a 

m 

> 


m 

A 

H 

m 

M 


> 


c  u  o 


B 

« 

O 

• 

4 

*4 

« 

X 

B 

9 

« 

o 

B 

e 

• 

4 

fN  IN 

5> 

e 

« 

X 

NN 

X 

4 

o 

X 

o 

• 

O 

XX 

M 

B 

• 

• 

•4 

•W 

u 

B 

o 

• 

9 

4  4 

o 

H 

• 

• 

4 

M  H 

u 

00 

« 

X 

m  m 

M 

B 

4 

« 

> 

o 

« 

B 

•  ^ 

NN 

B 

« 

O 

«  X 

B 

B 

fN  X 

• 

X 

«  B 

M  M 

9 

•  U 

• 

M 

4 

•  B 

• 

B 

4 

n  ■ 

B 

4 

^  ► 

• 

X 

•  M 

«  « 

M 

> 

4 

« 

O 

«  a 

M .«« 

•J 

0m  0m 

M.  ««> 

« 

M 

•  o  — 

M  M 

M  ^ 

4  4 

B 

X  9  O 

• 

X 

«  E  B 

w  w 

9  9 

M 

X 

>  4  • 

« 

M 

•  X  * 

M  «• 

n»  9 

O 

«l  B 

9 

X  o 

• 

O 

•  <*  o 

M  U 

o  o 

4 

>  w  . 

• 

X 

•  B 

«  • 

♦  -f 

B 

X  B  B 

O 

X  9 

B  O 

• 

o 

•  X  o 

M  A 

II 

BOO 

B 

IN  B  B 

O 

X 

• 

u 

«  M  IN 

X  X 

0^  0m 

> 

4  B  B 

B 

«  9  • 

O 

4 

• 

•  X  ** 

X  X 

B 

B  NX 

4 

•  B  -* 

X 

a 

e 

• 

X 

4 

•  H  X 

IN 

■  ■ 

■  ^  ^ 

«->  4 

x 

M 

• 

« 

M 

4 

«  9 

• 

o  o 

M  »• 

B 

*4  4 

B 

4  1'* 

M 

« 

o 

« 

x 

H 

«  M 

ij 

•J 

*4 

*4 

u  u 

9  9 

B 

0^  m0 

B 

'*  B 

M 

•N 

« 

M 

X 

•  X  X 

M 

M 

M 

1 

M 

M 

H 

•  MX 

X 

M  —  X 

B 

B 

1 

1 

1 

u 

« 

M 

«  H  O 

M 

n 

M 

M 

1  1 

119 

9 

4  9  9 

X  4  M 

9 

X 

M 

9 

« 

B 

X 

«  X  M 

• 

1 

I 

B 

B 

B 

•  • 

B 

'*  a 

U 

4 

0m 

X 

x 

• 

X 

H 

«  9  B 

lA 

w 

•  I 

^  ^  M 

M 

H 

^  0m  0m 

X 

•  B  X 

M 

X 

B 

B 

M 

a 

• 

M 

«  O  X 

fi4 

IN 

•  R 

X  X  H  H 

4 

B  4  4 

4 

X  ^ 

X 

X 

•4 

« 

X 

9 

« 

X 

B 

«  B  M 

M 

M  >• 

^  ■ 

■ 

U 

4 

9  B 

■ 

4 

H 

X 

X  o  • 

9 

M 

«  B  X 

R  R  X 

o 

1 

3  9 

M  X  O 

o 

9 

O  M  X 

9 

4  B  B 

O 

o  o 

w 

o 

B 

X  « 

o 

X 

•  9  M 

H 

H 

M 

O  M  o  O'  9  9  o 

v; 

U 

O  >  > 

U 

X  X  M  U 

0 

o 

o  4 

o 

B 

u  • 

B 

M 

«  B  O 

■J 

4 

B 

• 

B 

X 

IN  O 

4 

4 

o 

B 

B 

• 

9 

M 

«• 

• 

U 

V 

m 

B 

*N 

• 

• 

B 

O 

U 

u 

u  o 

W 

o  o 

u 

U  U 

u 

« 

u 

c 

c 

c 

u  o 

U 

OCriEC  lEITIAL  COaOITIOHS  rOR  SOLUTE  TRAESrORT 


DO  SS  N>J1.J2 
HT(N, I )>HeAO 
ii  COMTIHUE 


O 


Q 

O 

s 

■ 

»< 

•J 


h> 

< 

X 


vt 
—  D 

O  ^ 


« 

a 

—  O  O 

« 

• 

a 

o  o  o  >• 

a 

« 

a 

o  o  o 

o 

4 

• 

a 

O  aa  ^ 

o 

X 

• 

a 

«  ^  w  • 

o 

H 

« 

a 

w  Ik.  bl  • 

4a 

o 

« 

a 

a  u  U  ^  a 

w 

o  • 

« 

a 

a)  a  O  O  O 

01 

o  -* 

a 

a 

O  O  O  o 

O 

^  o 

a 

a 

M  o  ^ 

a 

•—» 

o 

a 

a 

>• 

—  • 

c 

o 

4 

a 

a 

O 

O  •  • 

cu 

o 

u  * 

a 

a 

O  —  O  M 

o 

a  o 

a 

a 

« 

• 

o  o  o 

» 

4  o 

a 

a  • 

X 

w 

^  o  o 

e 

a 

a  S 

u 

o 

W  o  o  »  • 

o 

fc* 

«.  r« 

a 

a 

•J 

a  ^  — 

o 

4 

— 

a 

M 

• 

W  o  o 

e 

X 

o  *4  • 

a 

a  M 

a 

H 

9  0  0  0 

a 

o  4  -« 

• 

a  O 

a 

*  w  u  o  o 

— « 

•j 

o  X  o 

a 

« 

—  Q  a  ^ 

a 

4 

^  H  o 

a 

a  a 

M 

* 

S  o  o  o - 

o 

— '  -4 

o 

a 

a  B 

B 

M 

o  a  a  ra  -4 

H 

» 

M  •  4 

o 

a 

a 

M 

o  9  >• 

a 

— «» 

U  -*  O 

o 

a 

»* 

* 

ra  • 

o 

e  o 

m 

a 

a  ^ 

a 

o 

u 

w  M  ^ 

• 

o 

«  o  o 

a 

a  »* 

a 

a 

>«  o  o  *  « 

o 

— .  -4  f4 

1*4 

a 

a 

a  M 

4  o 

o  o  — 

o 

4a 

O  4  •*-’ 

4 

•J 

M 

a 

a  a 

>• 

a 

u 

«  o  o  o  o 

o 

o  o  4 

— • 

a 

a 

o 

X 

—  ^  ^  o  o 

o 

o  o 

«• 

• 

• 

» 

a 

a 

o  '■>'**  o  o 

• 

M, 

^  o  • 

» 

• 

— • 

a 

a  M 

o  a  u  ^ 

a 

1*4 

<»<» 

O 

o 

o 

O 

e 

1*1 

a 

a  (H 

o  w  a  ' 

a 

a 

M  ^  O 

e 

o 

o 

o 

o 

o 

a 

a  a 

aa  Q  O  *4  -a 

o 

a  o  o 

o 

o 

e 

o 

o 

o 

o 

a 

M 

a  4 

^  O  O  9  M 

u 

4 

N  N,  o 

o 

M 

r4 

*a 

4«* 

a 

«9 

a  m 

a  a  a  \ 

s 

k*  M  a4 

r* 

«-* 

•» 

a 

a 

a  »< 

\  N  a 

X 

Q  01  w 

w 

H  *4 

«4 

H 

H 

a 

< 

a 

M  a  a 

2 

4  a  a 

a 

a 

n 

a 

9  9 

a 

m 

a  M 

Ha  4 

a 

u  4 

H 

a 

a  B 

u  o  a 

4 

H  a  a 

f«# 

a 

a  M 

a  a  s 

a 

01  H  O 

w 

a 

M 

a 

w  \ 

N 

\N  w 

H 

s 

a 

a 

a  3 

a  B  B 

a 

B  a  o> 

a 

a 

a 

H 

a  O 

O  O  0 

o 

O  ^  X 

X 

3 

a 

h 

a  m 

XX  X 

X 

n  X  a 

a 

h 

O  a 

D 

a  m 

XX  X 

X 

XXX 

c 

u 

B  « 

O 

a  3 

o  o  o 

o 

o  c  « 

tto 

a 

kl  a 

m 

a  M 

u  u  u 

u 

u  u  o 

a 

m 

a 

a 

a 

a 

a  a 

a 

a 

a 

a 

O 

a 

3 

a 

o 

a 

u 

M 

s 

I 

a 

B 

•< 

a 

u 


B 

M 

>« 


U  U  U  U  t)  o  o 


O 

VI 

a 

O 

H 

o 

•H 

>4 

M 

u 

H 

o 

M 

M 

H 

•4 

o 

a 

a 

O 

X 

o 

o 

M 

u 

4 

Q 

v4 

H 

H 

4 

M 

U 

u 

M 

hi  a 

u 

a 

a 

o 

a  ^ 

M 

H 

o 

•  a 

U 

>4 

a 

A 

a 

a 

a 

a  — 

M 

e 

4 

• 

^  o 

M 

e 

o 

«4 

a  o 

*4 

a 

1 

a  1-4  M 

a 

3 

a 

Q  *  3 

o 

lA 

o 

0  a 

a 

a 

o 

a  *-  M 

Q 

a4 

*-  O  H 

B 

4 

M 

4  a 

G 

M 

a 

o  w  u  o 

M 

a 

M 

Q 

M  a  u 

H 

a 

a 

M 

o 

0 

a 

u  u 

u 

u  li 

f 


V  u 


DEriHB  HEW  BOUHDASX  FLUXES 


ir  i Kooer ( N I . HE . 1 1  go  to  20 

LKB)  -  SQRT(  (XI  (M  )-X(N  )  )  • ‘2  ♦  (Yl(N)-Y(N)) 

L2(ll|  -  SQ«T( 1X2(N)-X(N) )**2  +  (Y2(N)-Y(N)) 
(2(H)  -  Q(N)  *  gl  (  N  )  'Ll  (  N  )/2  +  g2  (  N  )  ‘L2  (  N  )/2 
20  CONTINUE 


11  CONTINUE 


IW( 1 )«IEL 
IW( 2 )«JCL 
IW( 3 )»KEL 


o  ^  m  a 

O  MO  O 

^  H  H 

—  9  VI  VI 

M  M  ♦  ♦ 


kJ  I 

a—irto  lO 

>aw^^  m  ^  m 

—  •HH 

Ml  M  SMS 

h,  I  Nk.wOO*^0 
MnXMQOUOU 


-  < 

•  la 
<  X 
»* 

H  O  Ul  t 

Q  X  X  M 

O  ■<  H 

■  m 

■  XII 


I  I  n  ^  u 

M  *  D 

M  M  ■ 

O  O  -N^  M 

fn  m  a  «  H 

M  <  ■ 

O  O  S  X  o 

O  O  H  H  O 


M  H  M  »4 

O  w  0^1 

o  o 

^  ^  —  B  — 

»  ■  r*  9t  B  r> 

'  m  m  m 

t  M  ^  ^  M 

-I  ol  «u' 

»  M  ««  D  <«  M  •«  D  ■ 

^  B  I 

I  <«  W  M  u  M  M  . 

'  r*  H  H  H  m  H  5 

I  M  B  M  M  B  ( 

1  O  «  O  B  O  «  O  < 

:  O  X  U  X  O  X  U  I 


6  71  fOSMAT (/>  LErT-HASO-SXDE-KATBIX '/) 
719  POBMAT(/'BIOHt-HAMD-SlDE-MAIRIX’ /) 
917  FORMAT ( 4EI0 . 3 ) 


MULTIPLY  MATRIX  WITH  VECTOR  AT  THE  RIGHT-HAND  SIDE 


BOUIDARY  COHDZTIORS 


SUBEOOTIIIB  SBANO 


O' 


'  O  D 

'  o  m 


•j 

•j 

u 


u 

(*i 

Q 

O 

z 

o 

a 

V) 

z 


o 

a 

oJ 

8 

a 

M 

a 

o 

o 

a 

VI 

o 

#»1 

z 

M 

z 

• 

o 

o 

0* 

H 

8 

(J 

a 

Z 

M 

o 

o 

9 

k. 

a 

o 

a 

M 

VI 

w 

* 

u 

U 

8 

—  z 

z 

m 

O  8 

u 

a 

^  M 

a 

• 

V) 

O  — 

a 

« 

*  H 

•J 

a 

H 

8  a 

VI 

9  8 

a 

■J 

8  • 

a  « 

8 

o 

9 

•  — • 

u 

a 

*  -4 

8 

VI 

•4  >04 

a 

H 

1 

a 

IB 

u 

•  ♦  w  a  H 

9 

8  h> 

• 

D  D 

'  ■  8 

l/>  -MM 


H  H 
>88 

•  o  O 
I  u  u 

^  o 


z  z 

<  < 
u  u 

«  8 

\  \ 


a  a 


M  W 

a  u 


8  a 

o  o 
o  u 


ol 

o 


a 

VI 

8 


Z  Z  Z 

M  * - 


8  8 


ol  ^  U 
u  w  u 
•-»  a 


O  O 
>  o  o 


u  a ' 

a  M  w  ; 


a 
o  >— 

a  » 

M 

O  • 

Q  M 


Z  Z 

a  a 

o  a 

hi  (4 

•  IDS 
8  8 

M  M 

Z  Z  H  H 
—  —  8  8 
a  a  o  o 

9  9  t;  u 

o  m 


.  o  e 

O  o4  *4 

^  a  lo 
a  IN  «  M 


!^^^5 

(  a  a  M  H 
»  o  o  o  a  8 

»  o  o  o 
>000 


u  u  u  u  u  u 


w  o  u  u  u  o 


uuoouuouo 


*EAOS  IM  XVEBAGE  LIHEAB  VELOCITIES  FOB  COBCESTBATIOB  BUN 
IF  VELOCITIES  HEBE  BOT  CALCULATED  IB  SUBBOUTIHE  STEADY  OB  TBAHSI 


IP  ( J2.lt. MILXM) 


« 

A 

« 

A 

« 

o 

A 

• 

A 

• 

a 

0 

• 

A 

• 

4 

A 

« 

M 

• 

A 

« 

o 

•4 

• 

» 

a 

• 

O 

« 

o 

4 

« 

a 

e 

• 

u 

A 

• 

o 

o 

• 

A 

A 

« 

M 

o 

« 

• 

A 

o 

« 

O 

M 

« 

4 

• 

« 

M 

A 

• 

A 

o 

• 

9 

a 

• 

A 

• 

• 

«J 

»3 

• 

a 

• 

A 

•  ^ 

O 

• 

M 

•J 

• 

•  A 

A 

« 

u 

w 

« 

A 

«  iJ 

• 

a 

A 

• 

O 

•  4 

A 

A 

« 

o 

l« 

« 

A 

•  U 

O 

A 

• 

u 

« 

M 

• 

A 

U 

A 

« 

« 

A 

• 

A 

4  « 

B 

« 

A 

> 

« 

a 

«  M 

A 

O 

o 

• 

M 

»J 

• 

o 

•  O 

B 

U 

A 

« 

A 

w 

• 

M 

•  o  — 

O 

« 

4 

H 

D 

O 

« 

A 

«  a  o 

M 

« 

O 

« 

J 

> 

4 

• 

« 

M 

•  BO 

A 

f4 

o 

« 

9 

>  o 

• 

o 

•  o 

M 

•> 

•«« 

« 

U 

♦ 

» 

• 

a 

•  f« 

o 

« 

•4 

X 

H 

a> 

e 

• 

o 

•  u  — 

B 

M 

> 

• 

4 

M  rt 

4 

M 

«4 

m 

u 

«  M  A 

O 

O 

o 

• 

U 

■J  « 

O 

» 

A 

• 

•  m  mi 

u 

a 

r* 

o 

*4 

« 

m  « 

«• 

• 

o 

u 

«  M  4 

u 

a 

A 

• 

O 

a 

•J 

• 

« 

m 

4 

•  U 

ij 

A 

• 

a 

a 

» 

• 

a 

1 

•w 

m 

« 

a 

M 

•  u 

4 

o 

*4 

o 

» 

A 

« 

o 

A 

M 

• 

M 

A 

«  a  a 

M 

o 

u 

A 

M 

« 

o 

i  M 

M 

• 

*4 

•  «  o 

A 

M 

p># 

1 

M 

kJ 

P4 

« 

u  > 

.J 

M 

9 

• 

W 

a 

•  Am 

44 

9 

% 

X 

9 

« 

M 

A 

a 

A 

a 

a 

M 

«  9  A 

a 

a  A 

X 

a 

r* 

A  a 

« 

a 

O  1 

« 

> 

M 

4 

A  • 

M 

•  OB 

M 

M 

w 

A 

*4 

4  A 

« 

M 

« 

H 

«■» 

H 

A 

X 

9  • 

A 

A 

•  AM 

A 

O  A 

A 

A 

X  9 

« 

A 

9 

tf 

a 

t« 

A 

A  O  • 

9 

M 

•  AX 

M 

a 

4 

•4 

a 

A 

A  A 

o  • 

9 

O  4 

A 

k. 

o 

4 

o 

M  a  • 

O 

a 

•  9m 

a 

o 

M  O  4 

o 

A 

4 

O  U 

B  « 

O 

Q  > 

> 

M 

u 

K 

A 

A  A  « 

A 

M 

•  AO 

M 

u 

A  O 

u 

M 

X 

A  A 

M  • 

A 

A 

« 

A 

A 

A 

A 

« 

A 

e 

O 

e 

• 

9 

M 

M 

A 

A 

O 

O 

« 

9 

m 

A 

o 

• 

A 

O 

o 

A 

A 

A 

A 

• 

A 

•4 

• 

A 

14 

« 

u 

C 

c 

c 

c 

u  u  u 

c 

c 

c 

u  u 

V  u 

c 

c 

C 

C 

U 

O 

C 

C 

C 

u 

SUBftOUTZHE  COaCOCLeN,  VX,  VY,  VBCTB,  IBAND« 

A,  VNODE,  C,  CALT,  NIT,  NTM,  IT,  IW) 
COKMON/KODE/KODCU(  4  000  )  ,  KODCQ(  4000  ) ,  KODEP ( 4000  )  , 
KOD£C( 4  0  00  )  .  KOOCQ(  40  00  ) ,  KODCP  <  40  00  ) 


•  —  o  ■—  ftl 

•  O  ^  O  VI 

I  o  —  O  M 

>  O  V)  O  Q 

1^0  ^ 


—  ^  r.*  o  ‘ 

O  U  O  Oh  ' 
O  O  M 


■  •  U<J 

BOOO<<‘~'^  *«U  0 

M000S0<^0<<U  >  • 

oo'vfltaackteXi: 

»o 

—  H  — o 

orM>^0  '-'O'-'  oo 

OUUH-UI  Oi-*  0<N 

OOS-V)—  <J-  O*  O'-- 

r4  o  ‘-'•OH 

^  »  .  -o  -o—  — 

>«^^^o-->oo  -O  ►*  - 

000^0000  >  — 

-OOG—  '  O 

-'000*40'-  -OU  *0 

0^'»''«l»<<^-40r'l<  *-4 

O'—  —  '-S'-'^O— ’f*0O  * 

OW-«mO«MaaO*Jbl2oo 
r4UOO*J0Vir4i<SOOO 
-^Ci>4U<<\M'-XH<J«0 


M  \ 

N 

>•  Q  a 

H  N  H 

— ■  r4 

Q 

X 

N  O 

X 

Q  \ 

U 

H  — 

O 

< 

M  m 

4 

4  U 

X 

>  4 

■ 

0 

H  4 

4 

u  X 

M 

V) 

n 

M 

M  m 

4 

H  O 

|4 

Z 

M 

* 

« 

u  u 

«• 

v»  O 

u 

O 

M 

•4 

V  N. 

\  \ 

M 

M 

• 

1 

X  m 

■ 

X  X 

■ 

VI 

0 

M 

o  o 

O 

o  o 

o 

■ 

H 

X 

X  X 

X 

w 

4 

«« 

A 

X  2 

X 

X  X 

X 

X 

X 

o  o 

o 

o  o 

o 

M 

o 

o 

o  u 

u 

U  ^4 

u 

o 

O 

o 

o 

• 

• 

« 

«  • 

• 

0 

O  O  O  ^  H  • 
III  «  -  9  C 

—  ^  M  M  K 

n  n  n  M  o 

*  •  -  U  0  H  -4 

M  M  M  ■ 

- c  t  o  o 

4  O  0  »4  H  O  O 


4 

X 

X 

w 

X 

w 

o 

>• 

M 

u 

> 

0 

> 

• 

N 

« 

>• 

A 

M 

■«« 

»« 

X 

0 

X 

M 

> 

•V 

> 

H 

M 

4 

4 

VI 

X 

X 

X 

0 

0 

0 

bl 

*4 

0 

*4 

0 

4 

> 

4 

V) 

• 

M 

♦ 

— 

♦ 

a 

X 

M 

X  X 

X 

M 

X 

V) 

•»» 

«—• 

> 

w 

0 

VI  V) 

0 

• 

0 

M 

9  9 

I- 

H 

> 

0  0 

U 

<■« 

U 

VI 

0  0 

0 

X 

o 

X 

M  M 

>  ^ 

'V 

> 

4 

o  o 

N  X 

N 

N  X 

0 

'— 

4 

^4 

M 

♦  ♦ 

•  0 

X 

•  0 

*4 

«  9 

0 

•  9 

X 

O 

•  A 

—  0 

U 

->  0 

0 

■ 

X  X 

X  0 

4 

X  0 

4 

4 

•ip» 

M 

M 

X 

0  0 

M  O 

1 

0  O 

1 

01 

*4 

»•  H 

> 

> 

0 

4 

U  0 

•  ♦ 

•  ♦ 

■ 

0  0 

X 

A 

— 

o 

■ 

M 

>  > 

X  - 

K 

X  — 

X 

o 

«  « 

--  X  J 

-  X 

4 

9 

•4  '- 

4 

0 

*4 

■9 

H 

X  X 

4  0 

X 

9 

4  0 

% 

M 

>— •  ^ 

X  H 

0 

0 

X  0 

1 

U 

O 

•4  H 

0  0 

a 

0 

0  u 

9 

■ 

4  4 

•4  0 

4 

M 

«J  0 

mm. 

o 

X  X 

4  > 

'V 

o 

4  > 

* 

0 

*4 

0  0 

0 

1 

M 

*4  U 

I 

1 

1 

o 

O 

U 

4  4 

o 

♦ 

s 

— • 

— • 

•4 

0 

H 

•  1 

X 

X 

X 

0 

X 

u 

'V 

v> 

0 

9 

0 

•— 

0 

•4 

f>4 

M 

X 

X 

9 

H 

0  0 

«4 

•4 

1*4 

0 

M 

M 

mm 

4 

0  0 

0 

0 

0 

H 

0 

*4 

M  M 

0 

0 

0 

0 

X 

4 

♦ 

9 

o  a 

M 

M 

M 

O 

O 

9 

U 

O 

o 

o 

u 

4 

u 

*4 

« 

• 

• 

O 

4 

4 

X 

o 

> 

(9 

V 

0 

»4 

0 

•— 

J  *3  *4  *4  iJ  *4 

U  U  U  (tl  U  fal 

MM*')  «  M 


»4  U  *4  *4  *4  *4 

U  M  M  U  M  M 

M  M  M  M 


U  U  U  V  u  o 


KCL»KU(M) 


CVALUATE  AAEA  Of  ELCMEHT 


u  u 


( 

i 

I 


i 


« 

u  <«<* 

«  4 
^  u 

«  « 

<  N 


u 

n 

w 


u 


u 

I 


•J 

H 

n 

M 

U 

•  iJ  ^  U 

M  M  U 
—  M  afi 

X  I  «  • 

<<  *4  (N  m 


»  »  » 


14^ 
I  —  *J  ^ 

M  <J 

I  ^  » 

i  S  ^  >  M 

M  Ml 

'  f  O  I  M 

I  M  O  S 


n  u 

+  • 


—  X 


X  o 
— ’  ♦ 


o  m 
^  o 

♦  • 


4  ♦ 

I 


X 

M 


» 

9 


H 

V 

O 

s 


o 

■f 


u 

< 

H 


M 

X 

o 

o 

*■ 

J 

*4 

W 

y  M 

M  «• 

• 

ON  N 

t 

« 

X 

9  —  — 

-£  * 

M 

4 

M 

u 

TO 

M 

> 

MM  M 

«* 

« 

U 

M 

o 

N 

O 

9 

X  M  M 

M 

M  O  4 

M 

n 

M  O 

♦ 

M  ♦  ^ 

M 

O  B  M 

« 

O  B 

B 

On  n 

m 

y 

O  4  M 

BBS 

y 

4 

M 

O  4 

B  B 

y 

m 

* 

M  • 

< 

B  M  M 

X 

M 

B  M 

M 

•  M  M 

X 

-  •  I 

1 

M 

-  -  e 

z 

M 

O 

M 

M  .4  -» 

•  in 

M 

s 

n  M 

M 

m 

M  *4 

•  1  c;  ta 

M 

>« 

9 

•  ff4 

• 

M 

-•-43  ID 

X  ^  X  ^  ■ 

M  n  o  M  n  M 

M 

a4 

w  n  - 

M 

O  o  -» 

9 

M 

« 

-  a 

M  B 

M 

y 

i 

M  n  I  a 

—  B 

B  B  n  M 

m 

H 

M  1 

M  a4 
B  B 

M 

M 

m 

y 

-Ml*  -  M 

M  SMS 

M  O  O  ^  O 

M 

M 

r4  f>4  J 

0  o  S 

M 

M 

m  y 

4  B 

X  o 

i 

4 

*4  •  H 

M  B 

0  O  o 

> 

O 

M 

4 

B  M 

Q 

i 

M  M  9  U  M  U 

< 

o  Q  y 

4 

y  u 

O 

O  O  4  U 

B 

M  M 

— 

M 

M 

« 

Q 

1 

M  ^ 

M 

« 

o 

M 

m 

o 

§» 

M 

U 

U 

n 

M 

M 

o  u  u  u  u  o 


u  u  u 


TIMC-DCT 

■GO- ( MTIMEC-l )/MSTEP 


MULTlftt  HATRIZ  WITH  VECTOE  AT  THE  IICUT-HANO  SIDE 


Oh 

a 

o 

MR 

• 

»• 

a 

u 

o 

a 

4 

> 

a 

Q 

1 

a 

N 

u 

RM 

•H 

H 

a 

M 

t4 

4 

1 

M 

M 

U 

H 

a 

o  • 

M 

9 

> 

m 

a 

P4 

«• 

O  * 

♦ 

H 

a 

0 

»•  M 

U 

»« 

w 

«9 

9 

a 

4 

o  a 

w 

O 

o 

a 

«  o 

0  4 

a 

u 

1 

N  I* 

a 

M 

« 

o 

a 

MR 

—  p» 

Ri» 

u 

» 

4 

o 

a 

m4 

•MR 

M 

a 

♦ 

O  •f 

u 

u 

u 

X 

t» 

a 

o 

a 

a 

4 

M 

M 

a 

O 

a 

a 

o 

a 

MR 

a  — 

o 

a 

a 

a  o 

u 

ij 

n 

a 

a 

a 

a 

o 

a 

4 

a  M 

a 

a 

O 

•f 

o  a 

• 

•  • 

u 

9 

o 

4 

o 

o 

o 

4  a 

•  RM 

U 

Q 

o 

o 

0  4 

M 

« 

a 

« 

4 

•M 

a 

a 

a 

o 

a 

a 

H 

t*  ^ 

a 

a 

H 

a 

MR. 

M 

a  a 

R* 

o 

M 

o 

a 

M 

a 

a 

a 

M 

•RM 

wJ 

o  u 

a 

> 

a 

U 

a 

M 

M 

a  M 

i4 

h 

1 

« 

O 

o 

M 

R 

• 

o 

R 

o 

o 

* 

R 

1 

• 

•  a 

o 

1 

-9 

a 

R 

■RM 

a 

•  » 

4 

u 

*4 

4 

>• 

M 

p4 

• 

•4 

M 

■4 

M 

a 

a 

I* 

a 

•«  • 

> 

m4 

u 

H 

m4 

a 

u 

a  <j 

a 

• 

a 

1 

• 

o 

M 

1 

■ 

a 

#• 

• 

1 

♦ 

u 

u  •  a 

u 

1  o  a 

1 

■ 

oa 

4 

1  1 

H  a 

> 

H  a 

r* 

H 

H 

1 

9 

a 

9 

4 

M 

n 

n 

4 

4  9 

a 

a  . 

X 

9 

a 

M 

9 

a 

M  ^ 

i  9 

» 

N 

<wr 

4 

a 

MR 

a 

U 

• 

a 

a  a 

> 

a 

1 

a 

1  a 

1 

a 

« 

1 

1  1 

9 

M 

a 

a 

a 

n 

M 

a 

a 

M 

M 

o 

o 

a 

a 

a  M  M 

1 

M 

a 

H 

a 

O 

*•* 

4 

^  M 

O 

H 

1 

H 

m 

w 

H 

•J 

u 

^  R-  ** 

a 

•J 

H 

H 

ta 

->  H 

41 

-  H 

w 

o 

u  o 

M 

4 

M 

a 

*4 

a 

a 

4 

^  a 

4 

M 

a 

4 

M  a 

a 

M  a 

s 

o* 

o  0» 

w 

« 

o 

w 

o 

O 

u 

o 

o 

o 

8 

a 

a  u  o 

M 

o 

9 

O 

RM  o 

a 

O  O 

RM  O 

M 

»« 

» 

O 

Q 

O  4 

u 

a 

a 

u 

o 

o 

a 

M 

M  a  u 

4 

t< 

o  a 

u 

»4 

O 

a  o 

M 

a  Q 

4  U 

Iw 

■ 

1  1 

M 

a 

wJ 

»• 

M 

4 

•b 

M 

M 

M 

a 

a 

9 

o 

o 

a 

in 

> 

o 

a 

4k 

O 

a 

m  m 

» 

a 

pii 

X 

a 

•« 

M 

0 

a 

a 

o 

«n 

U 

u 

C 

C 

C 

C 

C 

c 

u 

u  u 

u 

w 

u  u 

u 

c 

c 

c 

u 

C 

c 

v 

u  u 

-A 


u  u  u 


BOUHDAHT  COHDITtOHS 


CALL  SUtftOUTZHS  SBAlfO 


o 

>0 

O  9 

H  a 

X 

o  •  — 

O  9  9 


<# 

r*  r* 


O 

O 

l« 

O 

O 

9 

9 

o 

H 

r» 

9 

9 

o 

O 

9  H 

9 

H 

■ 

H 

r«  9 

f>« 

9 

A 

» 

• 

m 

O 

e  M 

O 

9 

M 

O 

t-  9 

»« 

9 

9  »• 

O 

t* 

9 

9 

9  M 

9 

M 

9 

X 

M 

M 

9 

9 

♦ 

P 

H 

9 

•  • 

* 

—  O 

9 

«l 

M 

U  9 

U  9 

ift  -4 

M  a 

M 

9 

9  »•  — 

M  O 

»4  M 

O 

•  9 

9 

a  9 

O 

O 

■  9 

X  > 

*  X 

O 

• 

•1 

O  4 

O 

M  O 

A 

•«  *4  M  9 

hJ  t«  M 

9 

•4  • 

A 

H 

9  m 

9 

9 

9  l«  9 

U  ♦  t«  H 

9  9  t«  H 

in 

m  M 

B  9 

• 

9 

9  • 

O  t«  9  9 

M  O  9 

9 

• 

«  • 

O 

P 

9  •  *0 

•  >w 

U 

• 

D  — 

♦  9  •  • 

A  ♦  • 

o 

41  O 

9 

VI 

M  .14  • 

p4  9 

9  O  »« 

9  P  9  M 

9  9  9  9 

•4 

9  H  *4 

4 

9 

A 

1 

« 

o 

9 

■ 

•  9 

9  f 

M 

9  9  «4  X 

O  9 

o  o 

9 

9  X 

9 

o  o 

9 

A  9 

<4 

A 

A 

A 

hi 

»4 

1 

P 

9 

-*  a 

4  9 

P  P  t  «4 

9 

A 

P 

■  M 

. 

m 

P 

•» 

A 

4 

9 

9  9 

9 

9 

O  9  J  A 

1 

A  »• 

9 

4  A 

A 

A 

9 

A  A  A  A 

9 

A 

X 

A 

A 

M 

9 

M 

9  A 

O 

M  M  9  • 

A 

M 

M 

o 

M 

9  1 

M 

M 

O 

M 

4  4 

4 

4 

9 

a 

9 

1 

44 

• 

H 

H  H 

r* 

O 

«  A  9  9 

9 

w 

A 

A 

9  9 

<w 

A 

A 

X  X 

X 

X 

P 

•J 

9 

4 

M 

9 

mi  9 

X  9  9  4  X 

P 

9  4  X 

9 

A 

m  m 

9 

9 

A  O 

4 

9 

0 

0 

O 

o 

4  O 

A  O 

M 

^  O  9  M 

O 

A 

A 

O 

O 

9  *4 

A 

A 

O 

O 

4 

o  o 

O 

O 

9  B 

u 

9 

O 

o 

0  4  0 

o 

O  O 

14  O  A 

U  o  9  A 

9 

M 

M 

<9 

U 

9  A 

M 

M 

9 

u 

X 

A  A 

A 

A 

9  9 

4 

A 

9 

9 

9 

A 

A 

A 

e 

O 

O  O 

O 

e 

|4 

9 

m 

9 

A 

A 

A 

9  O 

O 

o 

SUtlOUTZVt  OUTrUT 


susAouTXMi  ourruT  (araAM,  aveL,  acoac,  acLca.  aaooE.  btime. 
a,  HT,  c,  vx,  V») 

oiMEasioa  a(2ooo),  iiT(2oao,<),  c(2oao,<i,  vxmooo),  vrMOOot 
coHHoa/PLOT/iPE ,  icoa,  cone,  coaiar,  xsize,  ysize.  xnax.  yka 


W*1TI  PLOTTiaa  PAXAMETEaS  foa  SOLUTE  TPAaspORT  OaLT 


SET  mOXTXVt  COaCEtfTftATXOVS  TO  SEKO 


30  COHTXMUE 


so  COMTIHUC 


■ACK  SUBSTITUTION 


so  coaTiauE 


tatatias44  LLLLLLLLLLLLLLt.LLLLLLLLLLLLLLLLL'.LLLLLLLLLLLLLLLLLLLLLLLtLlLLLLLLlLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLtL  4444444444 

aaaaaaaaaa  llllllllllllllllllllllll  ooiij^a  sma/xya  -  noT>»»oajoo  )a*adtnbi  tot&tq  llllllllllllllllllllllll  4444444444 

44444444  LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL  4444444444 


ril»  OUAl : iWCBNH.rtOaBAMlCOHOUC.DAT; 1  (977, (£,0),  last  raviasd  on  16-JUH-19t7  10:17,  is  a  15  block  saquantial  fils  OHBSd  by  tfZC 
( oaXBVtB ,WSIBBCBT I .  Tbs  cscords  srs  variabls  Isofth  with  isplisd  (CB)  csrtisgs  contcol.  Tbs  loa9sst  csco'rd  is  54  bytss. 


ao9C  (.ttt 


.  O  O 

oo 
o  o 

«  o 

. .  m  w4 

I  ««  lA  ^  m 

. .  •  •  o  o 

o  e 
o  o 


o  o 

e  o. 


«nor*^^'or<or>o<nooMr«p4 

#^r»oooooor*^oo^« 


•*o»^or*i«or‘w^r>‘0«^»^f4«< 

f'voo^r-oooo^or*^ 


lOOOOOOOOOOOOOOOO 


o 


o^^oo<nor»«»«Ao«4^r«r«»4 


0«l«l000«1int>^*4»4.4»4*4000 


oo  «»tf>«f«4X«lN9<#«a>O«00OO000 
oo  »^oo*^«nooooo<^otr»#4 


•  *  •  •  •  •  •  * 

.n*f^r*o-^OirtiAor*ooof^««o 

I 

^  t$  m . *  *oo  oooooMf*«^<N^'^r<f^<*< 

«or*^^inor»oo«oo<^ 
o  ^r«<>4tftoooor^«>%o<^<*« 


^OM00«4O*40«M0«^0<^e 


r>or»r^r»r>o<A.'^'N'VO*40o 

»*»ioooo«OMi^om«r‘r‘0 

oooo 


o 

Q  ^  . . 


•i«  o  m  Ml  o  r«  r« 
o  o  r«  r» 

o  «  iO  lA  lA  «>  w> 


000*0  000*0  1  t  T 

0  Lttx  X 

0  LtfT  »ty 

000*T  T  (If  «0« 


.000  0.000 


-40«>*0'40p4<0«««»«40««4i 


9 


«*ir<r>«4«i'«9999r«r4«n999ii*» 


•  9  -4  M  9  O  I 
>  *1  9  9  ^  9  « 


^'^^9*4999991 

90999r^r**i*-49« 

9r>f*•r*r•^>r^9991 


>r«94r  r>o<^9^f*>a, 

t999900000«-4. 


« 


OOOOOOA<9000<9<90 

oooooooooooooo 


o  o  o 
O  O  o 


I  lo*  l>«  ^  *4  ^  t 


>*40«40«40«40^0«^0-40«40»40«40«40««<9«40i>40«40«>«0«<*0«i««4a4 


O  O  O  O  O  O  o 
o  o  o  e  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 


oooooeeo 

oooooooo 

oooooooo 

oooooooo 

oooooooo 


'  o  o  o  o  o 


>000000 


•4*40<9*>^^»''•00••00«00(Nf*•OOOM>4000f«*f•00-•00►»0•4nP• 

rifN(tnooo«<'^or><»ooooo»<*4  9)^«iOf4«<oooo^or>r»*«««o^o  oooooooooooooo 

<*9«*«f>«<*<r4«4«<#'<^4<<vooooo<ir<<o^o«^^oooor*r»r^f>r>r»ooooo  f«r4«^*>^oo4niooor'r«oo 

f4*4#4»4p4<4r4<4«4<«4«4*4#4v4p«^44<4r4»4«4*4*4«4*4«4«4«4*4»4«4p4«4«4«4«4«4  4i4«4a4»4^«4«4«4«4p4>4i4w<«4 


eocoo'OT  ffot  tfor 
OODOO'Ol  CtOr  f9ST 
OOOOO'OT  9f€T  T96T 


I 

* 


I 


OfiHOOOOOOOO 

oooooooooo 

oooooooooo 

oooooooooo 

OOOOOOO'V'OO 


ooooooooeoo 

ooooooooooo 

oooaooooooo 

ooooooooooo 

eoooooooooo 


>00000000 


0*40*<0*40«40*40 


oooooo9oooooooiAr*^oo«>40o<*^r^ 

oo<«oooooooooM^ooooA«*^r«r«o«4o 

o  o 


o  oooo»-»or‘r»o-^ 


4 


